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ABSTRACT 

My project computed the one loop fermion self-energy for massless Dirac + 
Einstein in the presence of a locally de Sitter background. I employed dimen- 
sional regularization and obtain a fully renormalized result by absorbing all 
divergences with Bogliubov, Parasiuk, Hepp and Zimmermann (BPHZ) coun- 
terterms. An interesting technical aspect of my computation was the need 
for a noninvariant counterterm, owing to the breaking of de Sitter invariance 
by our gauge condition. I also solved the effective Dirac equation for mass- 
less fermions during inflation in the simplest gauge, including all one loop 
corrections from quantum gravity. At late times the result for a spatial plane 
wave behaves as if the classical solution were subjected to a time-dependent 
field strength renormalization of Z2{t) = 1 — ^GH"^ ln(a) + O(G^). I showed 
that this also follows from making the Hartree approximation, although the 
numerical coefficients differ. 
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1 INTRODUCTION 



My research focussed on infer how quantum gravity affects massless fermions 
at one loop order in the inflationary background geometry which corresponds 
to a locally de Sitter space. In the following sections, we will discuss what 
inflation is, why it enhances the effect of quantum gravity, how one can study 
this enhancement and why reliable conclusions can be reached in spite of the 
fact that a completely consistent theory of quantum gravity is not yet known. 

1.1 Inflation 

On the largest scales our universe is amazingly homogeneous and isotropic. 
It also seems to have nearly zero spatial curvature Based on these three 
features our universe can be described by the following geometry, 

ds"^ = —dt^ + a^{t)dx ■ dx . (1) 

The coordinate t is physical time. The function a{t) is called the scale factor. 
This is because it converts Euclidean coordinate distance y|| into physical 
distance a(t)||x — ^||. 

From the scale factor we form the redshift z{t), the Hubble parameter 
H{t) as well as the deceleration parameter q{t). Their definitions are: 

The Hubble parameter H{t) tells us the rate at which the universe is ex- 
panding. The deceleration parameter measures the fractional acceleration 
rate (a/a) in units of Hubble parameter. The current value of Hubble pa- 
rameter is , Hq = {7^tt)^^ — 2.3 X 10~^^Hz [IJ. From the observation of 
Type la supernovae one can infer go — —0.6 which is consistent with a 
universe which is currently about 30% matter and 70% vacuum energy. 

Inflation is defined as accelerated expansion, that is, q{t) < as well as 
H(t) > 0. During the epoch of primordial inflation the Hubble parameter 
may have been as large as Hj ~ 10^''Hz and the deceleration parameter is 
thought to have been infinitesimally greater than —1. The current values 
of the cosmological parameters are consistent with inflation, however, the 
phenomenological interest in my calculation concerns primordial inflation. 
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1.2 Uncertainty Principle during Inflation 

To understand quantum effects during inffation it is instructive to review the 
energy-time uncertainty principle, 

AEAt ~ 1 . (3) 

Consider the process of a pair of virtual particles emerging from the vacuum. 
This process can conserve 3-momentum if the particles have ±fc but it must 
violate energy conservation. If the particles have mass m then each of them 
has energy, 



E{k) = Vr«2+ II k P . (4) 

The energy-time uncertainty principle restricts how long a virtual pair of 
such particles with ±k can exist. If the pair was created at time t, it can last 
for a time At given by the inequality, 

2E{k)At ~ 1 . (5) 

The lifetime of the pair is therefore 

At = -^-^ . (6) 

— * 

One can see that in flat spacetime all particles with 7^ have a finite 
lifetime, and that massless particles live longer than massive particles with 
the same k. 

How does this change during inflation? Because the homogeneous and 
isotropic geometry shown by Equation [1] possesses spatial translation invari- 
ance it follows that particles are still labeled by constant wave numbers k, just 
as in fiat space. However, because k involves an inverse length, which must 
be multiplied by the scale factor a(t) to give the physical length, the physical 
wave number is k/a{t). Therefore the physical energy is not Equation H] but 
rather. 



E{t,k) = ^Jm^+\\k\\ya^{t) . (7) 
The left-hand side of the previous inequality becomes an integral: 

i-t+At 
Jt 



dt'2E{t', k) ~ 1 (8) 
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Obviously anything that reduces E{t', k) increases At. Therefore let us 
consider zero mass. Zero mass will simplify the integrand in Equation [S] to 
2|| A;||/a(t'). If the scale factor a{t) grows fast enough, the quantity 2||A;||/a(t') 
becomes so small that the integral will be dominated by the lower limit and 
the inequality of Equation [8] can remain satisfied even though At goes to 
infinity. Under these conditions with m = and a{t) = aje^^, Equation [S] 
gives, 

2 lU II ,^ ^HAt, < 1 _ 



Ha{t) ' 

From this discussion we conclude that massless virtual particles can live 
forever during inflation if they emerge with || A; ||~ Ha{t). 



1.3 Crucial Role of Conformal Invar iance 

One might think that the big obstacle to inflationary particle production is 
nonzero mass. However, the scale of primordial inflation is so high that a 
lot of particles are effectively massless and they nevertheless experience little 
inflationary production. The reason is that they possess a symmetry called 
"conformal invariance." 

A simple conformally invariant theory is electromagnetism in D = 4 
spacetime dimensions. Consider D dimensional electromagnetism, 

Cem = —F^pFp^g^PgP'^^g , (10) 

where F^^ = d^Ay — dyA^. Under a conformal transformation g'^^ = Vl?{x)g^i, 
and A'^ = the Lagrangian becomes, 

c = F^f^F.^n-^g'^pn-^g^'^n'^y^ = m""-^ (ii) 

Hence electromagnetism is conformally invariant m D = A. Other confor- 
mally invariant theories are the massless conformally coupled scalar, 

C = -^d,(j)d,<Pg^''V^ -^-{^^^R^ . (12) 

and massless fermions, 

C = i^e\j\zd,-^A,,dr'')^V^. (13) 
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Here 0' = Q^^~(f) and 'ip' = ^l^~tp under a conformal transformation. 

If the theory possesses conformal invariance, it is much more convenient to 
express the homogeneous and isotropic geometry of Equation [1] in conformal 
coordinates, 

dt = a(t)di] =^ ds^ = —dt^ + a^{t)dx ■ dx 

= o?{t){-d7f + dx ■ dx) . (14) 

Here t is physical time and rj is conformal time. In the (r^, x) coordinates, 
conformally invariant theories are locally identical to their flat space cousins. 
The rate at which virtual particles emerge from the vacuum per unit confor- 
mal time must be the same constant — call it F — as in flat space. Hence 
the rate of emergence per unit physical time is, 

dN _dNdT] _ T 
dt dr] dt a(t) 

One can see that the emergence rate in a locally de Sitter background is 
suppressed by a factor of 1/a (a ~ e^* ,H>0). Therefore any conformally 
invariant, massless virtual particles with || ^ ||~ Ha{t) can live forever but 
the problem is that they don't have much chance to emerge from the vacuum. 



1.4 Gravitons and Massless Minimally Coupled Scalars 

Not every massless particle is conformally invariant. Two exceptions are 
gravity and the massless minimally coupled (MMC) scalar, 

C = -^{R-2K)^g, (16) 

C = -U.cpd^g^^^g . (17) 

Here R is the Ricci scalar and A is the cosmo logical constant. From previous 
sections one can conclude that big quantum effects come from combining 

• Inflation; 

• Massless particles; and 

• The absence of invariance. 
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Therefore one can conclude that gravitons and MMC scalars have the po- 
tential to mediate vastly enhanced quantum effects during inflation because 
they are simultaneously massless and not conformally invariant. 

To see that the production of gravitons and MMC scalars is not sup- 
pressed during inflation note that each polarization and wave number behaves 
like a harmonic oscillator [5], H] , 

L = -mq — -mcu q , (18) 

with time dependent mass m(t) = a^{t) and frequency uj{t) = The 
Heisenberg equation of motion can be solved in terms of mode functions 
u{t, k) and canonically normalized raising and lowering operators and a, 

k^ 

q + 3Hq-\ — jQ' = =^ q{t) = u{t,k)a + u*{t,k)Q^ with [a,a^] = l, 

(19) 

The mode functions u{t^ k) are quite complicated for a general scale factor 
a{t) [5j but they take a simple form for de Sitter, 



H 



'2k^ 



ik 



u{t, k) = 1 - — — exp -— — . (20) 



Ha{t)i '^Ha{t) 



ik 



The (co-moving) energy operator for this system is, 

E{t) = \m{t)e{t) + \m{t)uj\t)q\t). (21) 

Owing to the time dependent mass and frequency, there are no stationary 
states for this system. At any given time the minimum eigenstate of E{t) 
has energy \uj{t), but which the state changes for each value of time. The 
state \VL) which is annihilated by a has minimum energy in the distant past. 
The expectation value of the energy operator in this state is, 

{v\E{t)\Q) = la^mt, k)\' + la{t)e\u{t, k)\\^ ^^^^^ = ^ + ^ • (22) 



If one thinks of each particle having energy k/a{t), it follows that the number 
of particles with any polarization and wave number k grows as the square of 
the inflationary scale factor, 

iV(a)^(^r- (23, 
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Quantum field theoretic effects are driven by essentially classical physics 
operating in response to the source of virtual particles implied by quantiza- 
tion. On the basis of Equation [23] one might expect inflation to dramati- 
cally enhance quantum effects from MMC scalars and gravitons, and explicit 
studies over a quarter century have confirmed this. The oldest results are 
of course the cosmological perturbations induced by scalar infiatons [B] and 
by gravitons [7]. More recently it was shown that the one- loop vacuum po- 
larization induced by a charged MMC scalar in de Sitter background causes 
super- horizon photons to behave like massive particles in some ways [H [U [10] . 
Another recent result is that the one-loop fermion self-energy induced by a 
MMC Yukawa scalar in de Sitter background reflects the generation of a 
nonzero fermion mass fHX \T2\ . 

1.5 Overview 

One naturally wonders how interactions with these quanta affect themselves 
and other particles. The first step in answering this question on the linearized 
level is to compute the one particle irreducible (IPI) 2-point function for the 
field whose behavior is in question. This has been done at one loop order 
for gravitons in pure quantum gravity [T3|, for photons [H [9] and charged 
scalars [H] in scalar quantum electrodynamics (SQED), for fermions pT] [T2] 
and Yukawa scalars [TB] in Yukawa theory, for fermions in Dirac + Einstein 
[16l| and, at two loop order, for scalars in 0^ theory [T7] . 

In the first part of my dissertation we compute and renormalize the one 
loop quantum gravitational corrections to the self-energy of massless fermions 
in a locally de Sitter background. The physical motivation for this exercise 
is to check for graviton analogues of the enhanced quantum effects seen in 
this background for interactions which involve one or more undifferentiated, 
massless, minimally coupled (MMC) scalars. Those effects are driven by the 
fact that inflation tends to rip virtual, long wavelength scalars out of the 
vacuum and thereby lengthens the time during which they can interact with 
themselves or other particles. Gravitons possess the same crucial property 
of masslessness without classical conformal invariance that is responsible for 
the inflationary production of MMC scalars. One might therefore expect a 
corresponding strengthening of quantum gravitational effects during infla- 
tion. 

Of particular interest to us is what happens when a MMC scalar is Yukawa 
coupled to a massless Dirac fermion for non-dynamical gravity. The one loop 
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fermion self-energy has been computed for this model and used to solve the 
quantum-corrected Dirac equation [TT], 



Powers of the inflationary scale factor a = e play a crucial role in under- 
standing this equation for the Yukawa model and also for what we expect 
from quantum gravity. The Yukawa result for the self-energy [Tlj consists 
of terms which were originally ultraviolet divergent and which end up, after 
renormalization, carrying the same number of scale factors as the classical 
term. Had the scalar been conformally coupled these would be the only con- 
tributions to the one loop self-energy. However, minimally coupled scalars 
also give contributions due to inflationary particle production. These are 
ultraviolet finite from the beginning and possesses an extra factor of aln(a) 
relative to the classical term. Higher loops can bring more factors of ln(a), 
but no more powers of a, so it is consistent to solve the equation with only 
the one loop corrections. The result is a drop in wave function which is 
consistent with the fermion developing a mass that grows as ln(a). A recent 
one loop computation of the Yukawa scalar self-mass- squared indicates that 
the scalar which catalyzes this process cannot develop a large enough mass 
quickly enough to inhibit the process [15]. 

Analogous graviton effects should be suppressed by the fact that the 
h^j^ipif) interaction of Dirac + Einstein carries a derivative, as opposed to the 
undifferentiated (pipip interaction of Yukawa theory. What we expect is that 
the corresponding quantum gravitational self-energy will consist of two terms. 
The most ultraviolet singular one will require higher derivative counterterms 
and will end up, after renormalization, possessing one less factor of a than the 
classical term. The less singular term due to inflationary particle production 
should require only lower derivative counterterms and will be enhanced from 
the classical term by a factor of ln(a). This would give a much weaker 
effect than the analogous term in the Yukawa model, but it would still be 
interesting. And note that any such effect from gravitons would be universal, 
independent of assumptions about the existence or couplings of unnaturally 
light scalars. 

The second part of my dissertation consists of using the IPI 2-point func- 
tion to correct the linearized equation of motion from Equation [2l] for the 
field in question. We employ the Schwinger-Keldysh formalism to solve for 
the loop corrected fermion mode function. In the late time limit we find that 




(24) 
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the one loop corrected, spatial plane mode functions behave as if the tree 
order mode functions were simply subject to a time-dependent field strength 
renormalization. The same result pertains for the Hartree approximation in 
which the expectation value of the quantum Dirac equation is taken in free 
graviton vacuum. 

1.6 The Issue of Nonrenormalizability 

Dirac + Einstein is not perturbatively renormalizable [18], however, ultravi- 
olet divergences can always be absorbed in the BPHZ sense [HI [201 ED [22] • 
A widespread misconception exists that no valid quantum predictions can 
be extracted from such an exercise. This is false: while nonrenormalizability 
does preclude being able to compute everything, that not the same thing as 
being able to compute nothing. The problem with a nonrenormalizable the- 
ory is that no physical principle fixes the finite parts of the escalating series 
of BPHZ counterterms needed to absorb ultraviolet divergences, order-by- 
order in perturbation theory. Hence any prediction of the theory that can be 
changed by adjusting the finite parts of these counterterms is essentially ar- 
bitrary. However, loops of massless particles make nonlocal contributions to 
the effective action that can never be affected by local counterterms. These 
nonlocal contributions typically dominate the infrared. Further, they cannot 
be affected by whatever modification of ultraviolet physics ultimately results 
in a completely consistent formalism. As long as the eventual fix introduces 
no new massless particles, and does not disturb the low energy couplings 
of the existing ones, the far infrared predictions of a BPHZ-renormalized 
quantum theory will agree with those of its fully consistent descendant. 

It is worthwhile to review the vast body of distinguished work that has 
exploited this fact. The oldest example is the solution of the infrared prob- 
lem in quantum electrodynamics by Bloch and Nordsieck [23], long before 
that theory's renormalizability was suspected. Weinberg [21] was able to 
achieve a similar resolution for quantum gravity with zero cosmological con- 
stant. The same principle was at work in the Fermi theory computation 
of the long range force due to loops of massless neutrinos by Feinberg and 
Sucher [25| [26] . Matter which is not supersymmetric generates nonrenormal- 
izable corrections to the graviton propagator at one loop, but this did not 
prevent the computation of photon, massless neutrino and massless, confor- 
mally coupled scalar loop corrections to the long range gravitational force 
[23 [2S1 [211 [SO]. More recently, Donoghue [SB [22] has touched off a minor 
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industry [531 EH ESI ESI EZ] by applying the principles of low energy effective 
field theory to compute graviton corrections to the long range gravitational 
force. Our analysis exploits the power of low energy effective field theory in 
the same way, differing from the previous examples only in the detail that 
our background geometry is locally de Sitter rather than fiat0 

2 FEYNMAN RULES 

When the geometry is Minkowski, we work in momentum space because 
of spacetime translation invariance. This symmetry is broken in de Sitter 
background so propagators and vertices are no longer simple in momentum 
space. Therefore we require Feynman rules in position space. We start from 
the general Dirac Lagrangian which is conformally invariant. We exploit 
this by conformally rescaling the fields to obtain simple expressions for the 
fermion propagator and the vertex operators. However, there are several 
subtleties for the graviton propagator. First of all, the Einstein theory is 
not conformally invariant. Secondly, there is a poorly understood obstacle to 
adding a de Sitter invariant gauge-fixing term to the action. We avoid this 
by adding a gauge-fixing term which breaks de Sitter invariance. That gives 
correct physics but it leads to the third problem, which is the possibility of 
noninvariant counterterms. Fortunately, only one of these occurs. 

2.1 Fermions in Quantum Gravity 

The coupling of gravity to particles with half integer spin is usually accom- 
plished by shifting the fundamental gravitational field variable from the met- 
ric g^v{x) to the vierbein e^m(x)|§ Greek letters stand for coordinate indices 
and Latin letters denote Lorentz indices, and both kinds of indices take val- 
ues in the set {0, 1, 2, ... , (D — 1)}. One recovers the metric by contracting 
two vierbeins into the Lorentz metric rf"^, 

Qf^uix) = e^,b{x)eyc{x)rf"' ■ (25) 

The coordinate index is raised and lowered with the metric {e^ b = g^^^vb)^ 
while the Lorentz index is raised and lowered with the Lorentz metric (e^ = 

"'^For another recent example in a nontrivial cosmology see D. Espriu, T. Multamaki 
and E. C. Vagenas, Phys. Lett. B628 (2005) 197, gr-q c/0503033, , 

^ For anot her approach see H. A. Weldon, Phys. Rev. D63 (2001) 104010, 
gr-qc/00090"86| 
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rf"^e^c)- We employ the usual metric-compatible and vierbein-compatible 
connections, 

gpa-i, = ^ = ^gP"" (^ga^,,u + gua,f, - gf,u,a) , (26) 

6/36;^ = =^ A^^d = e''^{e^d,t,-TP^^epd) . (27) 

Fermions also require gamma matrices, 7^^-. The anti-commutation rela- 
tions, 

{7^7'=}^(7V + 7V) = -2r/'"=/, (28) 

imply that only fully anti-symmetric products of gamma matrices are actually 
independent. The Dirac Lorentz representation matrices are such an anti- 
symmetric product, 

(29) 

They can be combined with the spin connection of Equation [27] to form the 
Dirac covariant derivative operator, 

V^^d^+ '-A^cdJ"" . (30) 
Other identities we shall often employ involve ant i- symmetric products, 

^b^c^d ^ ^[b^c^d] _ ^bc^d ^ ^db^c _ ^cd^b ^ ^^^^ 

I'r' = ^7NV^ + ^r/^V - V ■ (32) 

We shall also encounter cases in which one gamma matrix is contracted into 
another through some other combination of gamma matrices, 

7^, = -DI, (33) 

7V76 = (^-2)7% (34) 

^y/^^ = 47]"'^/- (D-4)7V , (35) 

7V7V7f> = 2YYl"+iD-A)YYY ■ (36) 

The Lagrangian of massless fermions is, 

-^^Dirac = i)e\-i^i'Df,i)^/^ . (37) 
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Because our locally de Sitter background is conformally flat it is useful to 
rescale the vierbein by an arbitrary function of spacetime a(x), 

e/3b = aepb =^ = a'^ e^^ . (38) 

Of course this implies a rescaled metric g^y, 

9,. = a^9^.u =^ g^'' = a-^r''- (39) 

The old connections can be expressed as follows in terms of the ones formed 
from the rescaled fields, 

rV = a-'{5^^a, + S^,a,-r''a^g,,)+T''^, (40) 
A^cd = -a~^{e^c'^^d-e^d^t^c)a,u + Af,cd ■ (41) 
We define rescaled fermion fields as follows, 

\1> = a^ip and = a^^p . (42) 



The utility of these definitions stems from the conformal invariance of the 
Dirac Lagrangian, 

/^Dirac = ^e'^b 7' , (43) 

where = d^ + ^A^cdJ"". 

One could follow early computations about flat space background [381 [39] 
in defining the graviton field as a first order perturbation of the (conformally 
rescaled) vierbein. However, so much of gravity involves the vierbein only 
through the metric that it is simpler to instead take the graviton field to be 
a first order perturbation of the conformally rescaled metric, 

dfiu = Vfiv + i^h^^ with = 16nG . (44) 

We then impose symmetric gauge {e^b = Cbp) to fix the local Lorentz gauge 
freedom, and solve for the vierbein in terms of the graviton, 

= (\/^^)/ = + \t^hpb - ^i^^h^h^b + • • ■ (45) 

It can be shown that the local Lorentz ghosts decouple in this gauge and one 
can treat the model, at least perturbatively, as if the fundamental variable 
were the metric and the only symmetry were diffeomorphism invariance 
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At this stage there is no more point in distinguishing between Latin letters for 
local Lorentz indices and Greek letters for vector indices. Other conventions 
are that graviton indices are raised and lowered with the Lorentz metric 
[h^^^ = rj^^hp^, h^^ = rj^fr]'^'^ hpcr) and that the trace of the graviton field is 
h = Tj^^hp^y. We also employ the usual Dirac "slash" notation, 



(46) 



It is straightforward to expand all familiar operators in powers of the 
graviton field, 
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b 



Z o 



(47) 
(48) 

-g = l + lnh + ln^h^ -]K^hf"'hp^ + ... (49) 
2 8 4 

Applying these identities to the conformally rescaled Dirac Lagrangian gives. 



^Dirac = ^^^^ + ^ hm^^-h^'^-f^ld^^-h^p^^^-f^^r'^^ 



-h'^--hP''h 

8 4 



pa 



4 8 P 



--hh 



HP, a 



~'~g^ p^vo-,^i ~^ fi^vp) ,a~^ a^^^P,l' 



^^A'jp-^ I + 0{k^) . (50) 



From the first term we see that the rescaled fermion propagator is the same 
as for flat space. 



iSj 



[x\ X 



r(f-i) 

Atx 2 



(51) 



where the coordinate interval is Ax^(x;x') = ||a; — x'p — {\Ti — ri'\ — i5Y . 

We now represent the various interaction terms in Equation EU] as vertex 
operators acting on the fields. At order k the interactions involve fields, 
and hap, which we number "1", "2" and "3", respectively. Each of 
the three interactions can be written as some combination Vj"^ of tensors, 
spinors and a derivative operator acting on these fields. For example, the 
first interaction is, 



(52) 
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Table 1: Vertex operators Uj^j'''^ contracted into jhaph^ 



# 


Vertex Operator 


# 


Vertex Operator 


1 




5 




2 




6 




3 




7 




4 




8 





Hence the 3-point vertex operators are, 

(53) 

The order interactions define 4-point vertex operators U^fj^'^ similarly, for 
example. 



The eight 4-point vertex operators are given in Table [H Note that we do not 
bother to symmetrize upon the identical graviton fields. 



2.2 The Graviton Propagator 

The gravitational Lagrangian of low energy effective field theory is, 

>CEinstein = T^{R " {D-2)k)^. (55) 

The symbols G and A stand for Newton's constant and the cosmological 
constant, respectively. The unfamiliar factor of D — 2 multiplying A makes 
the pure gravity field equations imply R^y = Kg^^ in any dimension. The 
symbol R stands for the Ricci scalar where our metric is spacelike and our 
curvature convention is, 

R = R^y = g^ {^''up,p ~ ^^pti,^ + ^^pa^ up ~ ^^va^ pfi)- (56) 

Unlike massless fermions, gravity is not conformally invariant. However, it 
is still useful to express it in terms of the rescaled metric of Equation [32] and 
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connection of Equation HUl 



-^Einstein = j^l^a""-^ R-2{D -l)a''-^r'' {^f^u -T" 

-(D-4)(D-l)a^-^^'^>^a^-(D-2)Aa^|y^ . (57) 

The factors of a which comphcate this expression are the ultimate reason 
there is interesting physics in this model! 

None of the fermionic Feynman rules depended upon the functional form 
of the scale factor a because the Dirac Lagrangian is conformally invariant. 
However, we shall need to fix a in order to work out the graviton propaga- 
tor from the Einstein Lagrangian in Equation [571 The unique, maximally 
symmetric solution for positive A is known as de Sitter space. In order to 
regard this as a paradigm for inflation we work on a portion of the full de Sit- 
ter manifold known as the open conformal coordinate patch. The invariant 
element for this is, 

ds"^ = o?(—drf + dx-dx\ where a{rf) = ——— , (58) 



and the /^-dimensional Hubble constant is if = ^A/{D — 1). Note that the 
conformal time 1] runs from —oo to zero. For this choice of scale factor we 
can extract a surface term from the invariant Lagrangian and write it in the 
form [¥T] . 

-^Einstein- Surface = {§-i)Ha''-'^r''r''hp„,,Ko + a^'-'^r^g'^r" 
^ ^\hap,^hpiyi, \hQ^pphfj^p -\- \hi^pph^piy i ho^pphpiyi, . (59) 

Gauge fixing is accomplished as usual by adding a gauge fixing term. 
However, it turns out not to be possible to employ a de Sitter invariant gauge 
for reasons that are not yet completely understood. One can add such a gauge 
fixing term and then use the well-known formalism of Allen and Jacobson 
[42j to solve for a fully de Sitter invariant propagator [l3l HHl HH HSl Wl\ . 
However, a curious thing happens when one uses the imaginary part of any 
such propagator to infer what ought to be the retarded Green's function of 
classical general relativity on a de Sitter background. The resulting Green's 
function gives a divergent response for a point mass which also fails to obey 
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the linearized invariant Einstein equation We stress that the various 

propagators really do solve the gauge-fixed, linearized equations with a point 
source. It is the physics which is wrong, not the math. There must be some 
obstacle to adding a de Sitter invariant gauge fixing term in gravity. 

The problem seems to be related to combining constraint equations with 
the causal structure of the de Sitter geometry. Before gauge fixing the con- 
straint equations are elliptic, and they typically generate a nonzero response 
throughout the de Sitter manifold, even in regions which are not future- 
related to the source. Imposing a de Sitter invariant gauge results in hy- 
perbolic equations for which the response is zero in any region that is not 
future-related to the source. This feature of gauge theories on de Sitter space 
was first noted by Penrose in 1963 |1H] and has since been studied for gravity 
[41j and electromagnetism [49] . 

One consequence of the causality obstacle is that no completely de Sitter 
invariant gauge field propagator can correctly describe even classical physics 
over the entire de Sitter manifold. The confusing point is the extent of the 
region over which the original, gauge invariant field equations are violated. 
For electromagnetism it turns out that a de Sitter invariant gauge can respect 
the gauge invariant equations on the submanifold which is future-directed 
from the source [50]. For gravity there seem to be violations of the Einstein 
equations everywhere [16]. The reason for this difference is not understood. 

Quantum corrections bring new problems when using de Sitter invariant 
gauges. The one loop scalar self-mass-squared has recently been computed 
in two different gauges for scalar quantum electrodynamics [13]. With each 
gauge the computation was made for charged scalars which are massless, 
minimally coupled and for charged scalars which are massless, conformally 
coupled. What goes wrong is clearest for the conformally coupled scalar, 
which should experience no large de Sitter enhancement over the flat space 
result on account of the conformal flatness of the de Sitter geometry. This 
is indeed the case when one employs the de Sitter breaking gauge that takes 
maximum account of the conformal invariance of electromagnetism in D = 
3 + 1 spacetime dimensions. However, when the computation was done in 
the de Sitter invariant analogue of Feynman gauge the result was on-shell 
singularities! Off shell one-particle-irreducible functions need not agree in 
different gauges [51] but they should agree on shell [52]. In view of its on- 
shell singularities the result in the de Sitter invariant gauge is clearly wrong. 

The nature of the problem may be the apparent inconsistency between de 
Sitter invariance and the manifold's linearization instability. Any propagator 
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gives the response (with a certain boundary condition) to a single point 
source. If the propagator is also de Sitter invariant then this response must be 
valid throughout the full de Sitter manifold. But the linearization instability 
precludes solving the invariant field equations for a single point source on the 
full manifold! This feature of the invariant theory is lost when a de Sitter 
invariant gauge fixing term is simply added to the action so it must be that 
the process of adding it was not legitimate. In striving to attain a propagator 
which is valid everywhere, one invariably obtains a propagator that is not 
valid anywhere! 

Although the pathology has not be identified as well as we should like, 
the procedure for dealing with it does seem to be clear. One can avoid the 
problem either by working on the full manifold with a noncovariant gauge 
condition that preserves the elliptic character of the constraint equations, 
or else by employing a covariant, but not de Sitter invariant gauge on an 
open submanifold [H]. We choose the latter course and employ the following 
analogue of the de Bonder gauge fixing term of flat space, 

Cgf = —a^'-Y'^F^F, , F^ = v'''{Kp,a - ^ + {D-2)Hah^X)- (60) 

Because our gauge condition breaks de Sitter invariance it will be neces- 
sary to contemplate noninvariant counterterms. It is therefore appropriate 
to digress at this point with a description of the various de Sitter symme- 
tries and their effect upon Equation EHl In our D-dimensional conformal 
coordinate system the ^D{D+1) de Sitter transformations take the following 
form: 

1. Spatial translations — comprising (D — l) transformations. 

r]' = 7] , (61) 
x'' = x' + e\ (62) 

2. Rotations — comprising ^{D — l){D — 2) transformations. 

T]' = T] , (63) 
= R^x^ . (64) 

3. Dilatation — comprising 1 transformation. 

T]' = krj , (65) 
x'' = kx\ (66) 
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4. Spatial special conformal transformations — comprising (D—l) trans- 
formations. 



77 = 



X 



V 



l-2e-x+\\9Px-x 



X* — 6'^x-x 



i-2e-x+\\ePx-x 



(67) 
(68) 



It is easy to check that our gauge condition respects all of these but the 
spatial special conformal transformations. We will see that the other sym- 
metries impose important restrictions upon the BPHZ counterterms which 
are allowed. 

It is now time to solve for the graviton propagator. Because its space and 
time components are treated differently in our coordinate system and gauge 
it is useful to have an expression for the purely spatial parts of the Lorentz 
metric and the Kronecker delta, 



V 



ViJ^u + Sl6l and 



(69) 



The quadratic part of £Emstcin + -^^gf can be partially integrated to take the 
form ^h'^^D^J'^hpcr, where the kinetic operator is, 



— I 2 



[^^"^ 2(^-3)'''^°° 



and the three scalar differential operators are. 



Da 

Dc , 



-V D\D- 



(70) 

(71) 
(72) 

(73) 



The graviton propagator in this gauge takes the form of a sum of constant 
index factors times scalar propagators. 



(74) 
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The three scalar propagators invert the various scalar kinetic operators, 

Dj xiAi{x;x') =i6^{x-x') for I = A,B,C, (75) 
and we will presently give explicit expressions for them. The index factors 



are, 



rjiA 

tif-'- pa 

rpB 

fJ'i^-'- pa 

rpC 

p.i^''' pa 



Ipv'lpa 5 



'{p'h){p^a) , 

2 



(76) 
(77) 



(D-2)iD-3) 



(D-3)5X + Vpu iD-3)5X + Vpa -(78) 



With these definitions and Equation [75] for the scalar propagators it is straight- 
forward to verify that the graviton propagator of Equation [71] indeed inverts 
the gauge-fixed kinetic operator, 



D 



pa 



pv 



X I 



pa^ 



5^;^5^\5''{x-x' 



(79) 



The scalar propagators can be expressed in terms of the following function 
of the invariant length ^{x\ x') between x^ and x'^^ 



aa'H^{\ 



X — X 



■/||2 



= aa' H'^ Ax'^ {x\ x') 
(|?7 — r/'l — z5)^ 



(80) 
(81) 



The most singular term for each case is the propagator for a massless, con 
formally coupled scalar |53j . 



iAcfix; x') 



(4vr; 



V 2 



1 



4xf-i 



(82) 



The A-type propagator obeys the same equation as that of a massless, min- 
imally coupled scalar. It has long been known that no de Sitter invariant 
solution exists [M] • If one elects to break de Sitter invariance while preserv- 
ing homogeneity of Equations [6T1I62] and isotropy of Equations [631IM] — this 
is known as the "E(3)" vacuum [55j — the minimal solution is |56l 151 

iAA{x; x) = iAci{x; x) 



T{D-1) J D r 



2 ' 



(47r)f r(f) \D-AT{D-r 

W^kh r(n+f) 



^ ^ — TT cot (— -D) + In(aa') 



2 

T{n + § + l] 



n-f + 2 r(n + 2) 



^+2 



.(83) 
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Note that this solution breaks dilatation invariance of Equations EMM] in ad- 
dition to the spatial special conformal invariance of Equations 1^711^ broken 
by the gauge condition. By convoluting naive de Sitter transformations with 
the compensating diffeomorphisms necessary to restore our gauge condition 
of Equation [60] one can show that the breaking of dilatation invariance is 
physical whereas the apparent breaking of spatial special conformal invari- 
ance is a gauge artifact |5$j . 

The B-type and C-type propagators possess de Sitter invariant (and also 
unique) solutions, 



iAB{x;x') =zArf(x;x') - ^— irE 1 da^^ (t)" 

(47r)-^ol r(n + y) 



tAc{x;x')=tA,i{x;x') + —-^Yl (^ + 1) , d 



n 



(4vr)f„t^,r r(n+f) V4. 

-in + 3) , ^ , M -) ' \. 85 

V 2 ' r(n+2) J 

They can be more compactly, but less usefully, expressed as hypergeometric 
functions [591 160] . 



(47r) 2 i (yj I ^ 



. , ,^ rm-3)r(2) / d 



(47r)T El 



These expressions might seem daunting but they are actually simple to use 
because the infinite sums vanish in D = 4, and each term in these sums 
goes like a positive power of y[x] x'). This means the infinite sums can only 
contribute when multiplied by a divergent term, and even then only a small 
number of terms can contribute. Note also that the -B-type and C-type 
propagators agree with the conformal propagator in D = 4. 

In view of the subtle problems associated with the graviton propagator 
in what seemed to be perfectly valid, de Sitter invariant gauges [161[1T], it is 
well to review the extensive checks that have been made on the consistency 
of this noninvariant propagator. On the classical level it has been checked 
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that the response to a point mass is in perfect agreement with the hnearized, 
de Sitter-Schwarzchild geometry [H]. The hnearized diffeomorphisms which 
enforce the gauge condition have also been exphcitly constructed [61]. Al- 
though a tractable, D-dimensional form for the various scalar propagators 
iAi{x]x') was not originally known, some simple identities obeyed by the 
mode functions in their Fourier expansions sufficed to verify the tree order 
Ward identity [^I] • The full, D-dimensional formalism has been used recently 
to compute the graviton 1-point function at one loop order [62]. The result 
seems to be in qualitative agreement with canonical computations in other 
gauges [631 El]- A D = 3 + l version of the formalism — with regularization 
accomplished by keeping the parameter 5 7^ in the de Sitter length function 
y{x; x') Equation [HD — was used to evaluate the leading late time correction 
to the 2- loop 1-point function [65l [66]. The same technique was used to 
compute the unrenormalized graviton self-energy at one loop order [13]. An 
explicit check was made that the flat space limit of this quantity agrees with 
Capper's result [67] for the graviton self-energy in the same gauge. The one 
loop Ward identity was also checked in de Sitter background [13]. Finally, 
the D = A formalism was used to compute the two loop contribution from 
a massless, minimally coupled scalar to the 1-graviton function [68] • The 
result was shown to obey an important bound imposed by global conformal 
invariance on the maximum possible late time effect. 

2.3 Renormalization and Counterterms 

It remains to deal with the local counterterms we must add, order-by-order in 
perturbation theory, to absorb divergences in the sense of BPHZ renormaliza- 
tion. The particular counterterms which renormalize the fermion self-energy 
must obviously involve a single ip and a sing le At one loop order the 
superficial degree of divergence of quantum gravitational contributions to 
the fermion self-energy is three, so the necessary counterterms can involve 
zero, one, two or three derivatives. These derivatives can either act upon 
the fermi fields or upon the metric, in which case they must be organized 
into curvatures or derivatives of curvatures. We will first exhaust the pos- 
sible invariant counterterms for a general renormalized fermion mass and a 

•^Although the Dirac Lagrangian is conformally invariant, the counterterms required to 
renormaUzc the fermion self-energy wiU not possess this symmetry because quantum grav- 
ity does not. We must therefore work with the original fields rather than the conformally 
rescaled ones. 
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general background geometry, and then specialize to the case of zero mass 
in de Sitter background. We close with a discussion of possible noninvariant 
counterterms. 

All one loop corrections from quantum gravity must carry a factor of ~ 
mass~^. There will be additional dimensions associated with derivatives and 
with the various fields, and the balance must be struck using the renormalized 
fermion mass, m. Hence the only invariant counterterm with no derivatives 
has the form, 

K^m^^fjipy/^ . (88) 

With one derivative we can always partially integrate to act upon the ip field, 
so the only invariant counterterm is, 

K'^m^i)ipi)^f^ . (89) 

Two derivatives can either act upon the fermions or else on the metric to 
produce curvatures. We can organize the various possibilities as follows, 

^^^m'^lJ{^]^Y'^l)^J—g , K^mRipipyZ—g . (90) 

Three derivatives can be all acted on the fermions, or one on the fermions 
and two in the form of curvatures, or there can be a differentiated curvature, 

n\^^R^^^-^g^'^R)^^"^^V,^P^ , K^e'^^i^^T^^v^ • (91) 

Because mass is multiplicatively renormalized in dimensional regulariza- 
tion, and because we are dealing with zero mass fermions, counterterms in 
Equations [HHl [89] and [90] are all unnecessary for our calculation. Although all 
four counterterms of Equation [HI] are nonzero and distinct for a general met- 
ric background, they only affect our fermion self-energy for the special case 
of de Sitter background. For that case Rf^i, = {D — l)H^g^u, so the last two 
counterterms vanish. The specialization of the invariant counter-Lagrangian 
we require to de Sitter background is therefore, 

R 

A£i„v = aiK^i^[{ti;f + ^^^—^)ti;^y^ + a2K^R^iVi^^ , (92) 
— > aiK,^^{i^a~H^a~^+ lO ^^^ + a2{D-l)DK,^H^m^^ . (93) 
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Here ai and 02 are D-dependent constants which are dimensionless for D = A. 
The associated vertex operators are, 



Of course Ci is the higher derivative counterterm mentioned in section 1. It 
will renormalize the most singular terms — coming from the iAd part of the 
graviton propagator — which are unimportant because they are suppressed 
by powers of the scale factor. The other vertex operator, C2, is a sort of 
dimensionful field strength renormalization in de Sitter background. It will 
renormalize the less singular contributions which derive physically from in- 
flationary particle production. 

The one loop fermion self-energy would require no additional countert- 
erms had it been possible to use the background field technique in background 
field gauge [^11 EOl ED E2]- However, the obstacle to using a de Sitter in- 
variant gauge obviously precludes this. We must therefore come to terms 
with the possibility that divergences may arise which require noninvariant 
counterterms. What form can these counterterms take? Applying the BPHZ 
theorem |19l |20l [21], [22] to the gauge-fixed theory in de Sitter background 
implies that the relevant counterterms must still consist of times a spinor 
differential operator with the dimension of mass-cubed, involving no more 
than three derivatives and acting between \1/ and As the only dimension- 
ful constant in our problem, powers of H must be used to make up whatever 
dimensions are not supplied by derivatives. 

Because dimensional regularization respects diffeomorphism invariance, it 
is only the gauge fixing term in Equation [001 that permits noninvariant coun- 
terterms]^ Conversely, noninvariant counterterms must respect the residual 
symmetries of the gauge condition. Homogeneity of Equations [6T1I62I implies 

^One might think that the they could come as well from the fact that the vacuum 
breaks de Sitter invariance, but symmetries broken by the vacuum do not introduce new 
counterterms [73]. Highly relevant, explicit examples are provided by recent computations 
for a massless, minimally coupled scalar with a quartic self-interaction in the same locally 
de Sitter background used here. The vacuum in this theory also breaks de Sitter invariance 
but noninvariant counterterms fail to arise even at two loop order in either the expectation 
value of the stress tensor [SH] [ST] or the self- mass-squared [T7] . It is also relevant that 
the one loop vacuum polarization from (massless, minimally coupled) scalar quantum 
electrodynamics is free of noninvariant counterterms in the same background [9]. 
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that the spinor differential operator cannot depend upon the spatial coor- 
dinate X*. Similarly, isotropy of Equations I^HIIMI requires that any spatial 
derivative operators di must either be contracted into 7* or another spatial 
derivative. Owing to the identity, 



{I'd, 



(96) 



we can think of all spatial derivatives as contracted into 7*. Although the 
temporal derivative is not required to be multiplied by 7° we lose nothing by 
doing so provided additional dependence upon 7° is allowed. 

The final residual symmetry is dilatation invariance shown by Equations 
I^M^ It has the crucial consequence that derivative operators can only 
appear in the form a~^d^. In addition the entire counterterm must have an 
overall factor of a, and there can be no other dependence upon 77. So the most 
general counterterm consistent with our gauge condition takes the form. 



(97) 



where the spinor function S{b, c) is at most a third order polynomial function 
of its arguments, and it may involve 7° in an arbitrary way. 

Three more principles constrain noninvariant counterterms. The first of 
these principles is that the fermion self-energy involves only odd powers of 
gamma matrices. This follows from the masslessness of our fermion and 
the consequent fact that the fermion propagator and each interaction vertex 
involves only odd numbers of gamma matrices. This principle fixes the de- 
pendence upon 7° and allows us to express the spinor differential operator 
in terms of just ten constants Pi, 



f^H'^aS[{Ha)-^-f°do,iHa)-Wdi) = fi:'a|/?i(a- V^o)^ 

+/32[(a-S%)'(a-V5i)] + /33[ia-Vdo)ia-Wdrf] + (3,ia-W^^f 

+i/2||/?8(a-S%) +/59(a-V5.)) +^V/3io|. (98) 

In this expansion, but for the rest of this section only, we define noncommut- 
ing factors within square brackets to be symmetrically ordered, for example, 
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+ l{a-Vdo){a-'7%){a-'j'do) + ^(a" V5^)(«"'7%)' • (99) 

The second principle is that our gauge condition of Equation EDI becomes 
Poincare invariant in the flat space limit of if — > 0, where the conformal time 
is ?7 = —e~^^/H with t held fixed. In that limit only the four cubic terms of 
Equation [98] survive, 



(100) 



Because the entire theory is Poincare invariant in that limit, these four terms 
must sum to a term proportional to (7^(9^)^, which implies. 



Pi = = = Pa . 



(101) 



But in that case the four cubic terms sum to give a linear combination of the 
invariant counterterms of Equation [91] and Equation [ 



K'a{ {a-'i^doY + 3 (a-^-f^doYia-^Ydi) 



+3 



(a-S%)(«"'7''9i)' + {a~'Y^^Y } = n'^a-'^a-'^ . (102) 



Because we have already counted this combination among the invariant coun- 
terterms it need not be included in S. 

The final simplifying principle is that the fermion self-energy is odd under 
interchange of and x'^, 



— I 



[x: X 



+i 



[X :x) 



(103) 



This symmetry is trivial at tree order, but not easy to show generally. More- 
over, it isn't a property of individual terms, many of which violate Equation 
11031 However, when everything is summed up the result must obey Equation 
11031 hence so too must the counterterms. This has the immediate conse- 
quence of eliminating the counterterms with an even number of derivatives: 
those proportional to P^^-r and to Piq. We have already dispensed with /3i_4, 
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which leaves only the linear terms, /?8-9- Because one linear combination of 
these already appears in the invariant of Equation [HS] the sole noninvariant 
counterterm we require is, 

ACr^on = "^Ca* where Csij = asK^HH'^ij . (104) 



3 COMPUTATIONAL RESULTS 

For one-loop order the big simplification of working in position space is that 
it doesn't involve any integrations after all the delta functions are used. 
However, even though calculating the one loop fermion self energy is only 
a multiplication of propagators, vertices and derivatives, the computation is 
still a tedious work owing to the great number of vertices and the complicated 
graviton propagator. Generally speaking, we first contract 4-point and pairs 
of 3-point vertices into the full graviton propagator. Then we break up the 
graviton propagator into its conformal part plus the residuals proportional 
to each of three index factors. The next step is to act the derivatives and 
sum up the results. At each step we also tabulate the results in order to 
clearly see the potential tendencies such as cancelations among these terms. 
Finally, we must remember that the fermion self energy will be used inside 
an integral in the quantum-corrected Dirac equation. For this purpose, we 
extract the derivatives with respect to the coordinates "x'^" by partially 
integrating them out. This procedure also can be implemented so as to 
segregate the divergence to a delta function that can be absorbed by the 
counterterms which we found in chapter 2. 



3.1 Contributions from the 4-Point Vertices 

In this section we evaluate the contributions from 4-point vertex operators of 
Table [H The generic diagram topology is depicted in Figure [H The analytic 
form is, 



4pt 



7=1 



') = H^^7ii ^ alB^pa [x] x') 6^ {x - x') . 



(105) 



And the generic contraction for each of the vertex operators in Table [T] is 
given in Table [2J 
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Figure 1: Contribution from 4-point vertices. 
Table 2: Generic 4-point contractions 



I 


iU\,]{x-x')iUf'''5^{x-x') 


1 


-\K^i[\APp]{x]x) ^S^{x-x') 


2 


i/«2z[-/3A„^](a;;a;) ^5^{x-x') 


3 


i/t^ z[°^Ap^](a;; x) Yd" 


4 


-f/t^ i5A„^](x; x) 7^9-^ 


5 


d'^iWd.p^]{x; x') YJ"^ 5^{x-x') 


6 


9;2["^A„,](x; x') YJ^" 5^(x-a;') 


7 


d^i[''pA^„]{x- x) -f^J''^' 5'^{x-x') 


8 


in^ d'(^i[afsAp^]{x; x') YJ'"" S^{x~x') 



From an examination of the generic contractions in Table [2] it is apparent 
that we must work out how the three index factors [a/jT^o-] which make up the 
graviton propagator contract into t]"'^ and t]"''. For the A-type and -B-type 
index factors the various contractions give, 



rpA 



D-3 

rpB' 

aP''- pa 



Vpa 
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^ap 



^ap 



rpA 

o^pJ-pa 
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(106) 



28 



For the C-type index factor they are, 



al3 



T. 



V 



ap 



fjiC 



' {D-2){D-3) ' 



(108) 



{D-2){D-3) 

On occasion we also require double contractions. For the A-type index factor 
these are, 



rpA 

aP-l-pa 



afi-'- pa 



D{D-l)-2 



D-1 
D-3 



(109) 



The double contractions of the -B-type and C-type index factors are, 



^a/S^pa 



af3 
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pa 
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all''- pa 



^ap^l3a 
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afSJ-pa 



rjiC 



2{D-l) , (110) 
(D2-5D + 8) 



\D~2){D-?>) 



.(111) 



(D-2)(D-3) 

Table [3] was generated from Table [2] by expanding the graviton propagator 
in terms of index factors, 



a/B^pa 



(x: X 



rpA 

a/3-' pa 



zA^(x; x') + 



rpB 

a/3-' pa 



iAsix: x) + 



Tc 



pa 



iAc(x:x') 



(112) 

We then perform the relevant contractions using the previous identities. Re- 
lation [32] was also exploited to simplify the gamma matrix structure. 

From Table [3] it is apparent that we require the coincidence limits of zero 
or one derivatives acting on each of the scalar propagators. For the A-type 
propagator these are, 

HD-2 y(D - 1) 



lim iA^(x;x') 

x'— >x 

lim dfj,iAA{x; x') 



(47r)f r(f) 
(47r)f r(f) 
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(114) 



The analogous coincidence limits for the S-type propagator are actually finite 
in D = 4 dimensions. 



lim iAsix; x') 
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Table 3: Four-point contribution from each part of the graviton propagator. 
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Figure 2: Contribution from two 3-point vertices. 



The same is true for the coincidence hmits of the C-type propagator, 
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Our final resuh for the 4-point contributions is given in Table HI It was 
obtained from Table [3] by using the previous coincidence limits. We have also 
always chosen to re-express conformal time derivatives thusly, 



;ii9) 



A final point concerns the fact that the terms in the final column of Table H] 
do not obey the reflection symmetry. In the next section we will find the 
terms which exactly cancel these. 



3.2 Contributions from the 3-Point Vertices 

In this section we evaluate the contributions from two 3-point vertex oper- 
ators. The generic diagram topology is depicted in Figure [21 The analytic 
form is, 



3pt 
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Table 4: Final 4-point contributions. All contributions are multiplied by 
^2gD-2r(D-i) ^ ^gg^g ^ ^ |cot(^)-ln(a). 
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Table 5: Generic contributions from the 3-point vertices. 



I 


J 


iV7'^(x) i[S]{x; x') iV^^ix') i[„/3Ap^](x; x') 


1 


1 




[X]X) 


1 


2 


-lK'^Yd''6'^{x-x') z["„Ap^](x; a;) 


1 


3 


liK'^YJ''''^^ix-x')d'^i[' 


aApo- (x, X ) 


2 


1 


jK^d'^A7'^d'^i[S]{x; x')Yi 
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2 
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~\t^d'P{YdH[S]{x-x')-i'' 
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2 


3 


-\iK'^ -i^dl^i[S] (x; x') YJ'"'"d'^ i 
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3 


1 


-\iK'^dl{YJ^^ i[S]{x; x') Yd^ i 


U^pp]{x-x')} 


3 


2 


\tK'^d'P{YJ'^''i[S] {x; x') Ydf, i[c 


yf3^p„]{x] X')} 


3 


3 


-\K^YJ^^i[S]{x]x')YJ"''^^.^'y 


i[ap/^pa\{x; X') 



Because there are three 3-point vertex operators of Equation [531 there 
are nine vertex products in Equation 11201 We label each contribution by the 
numbers on its vertex pair, for example, 



I-J 



= X i S {x;x') X iVf{x') X i 



[X\ X 



'1211 



Table [5] gives the generic reductions, before decomposing the graviton prop- 
agator. Most of these reductions are straightforward but two subtleties de- 
serve mention. First, the Dirac slash of the fermion propagator gives a delta 
function. 



s 



\X\ X 



i6 (x — x') . 



This occurs whenever the first vertex is 1 = 1, for example. 



1-3 



■r]^^i(^ X i 



S 



IK 



(x; x') X Yr^d'^ X I 



(x: X 



IK 



Yr^5'^{x-x')d'i 



^ A 



x; X 



(122) 

(123) 
(124) 



The second subtlety is that derivatives on external lines must be partially 
integrated back on the entire diagram. This happens whenever the second 
vertex is J=l or J = 2, for example. 



2-2 



2 ' 



S 



IK 



(x;x') X -—Yid'^^t X i 



afi^pa 



(125) 
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s 



[x; x) 



..(126) 



In comparing Table [5] and Table [2] it will be seen that the 3-point contri- 
butions with 1=1 are closely related to three of the 4-point contributions. 
In fact the [1—1] contribution is —2 times the 4-point contribution with / = 1; 
while [1—2] and [1—3] cancel the 4-point contributions with / = 3 and / = 5, 
respectively. Because of this it is convenient to add the 3-point contributions 
with / = 1 to the 4-point contributions from Table HI 
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+ 
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4V^ D-3 
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D- 



D-3 



aif7°U^(a;-a;').(127) 



In what follows we will focus on the 3-point contributions with / 
/ = 3. 



2 and 



3.3 Conformal Contributions 

The key to achieving a tractable reduction of the diagrams of Fig. 2 is that 
the first term of each of the scalar propagators iAj{x]x') is the conformal 
propagator iAcf{x; x). The sum of the three index factors also gives a simple 
tensor, so it is very efficient to write the graviton propagator in the form, 



fiuApo- 



[x: X 



D-2 



T]p^T]p„'^iAci{x] x') 

+ J2 [t^uT^i6AI{x■,x 



I=A,B,C 



;i28) 



where i6Aj{x; x') = iAj{x; x')—iAcf{x; x'). In this subsection we evaluate the 
contribution to Equation 11201 using the 3-point vertex operators of Equation 
53] and the fermion propagator of Equation [51] but only the conformal part 
of the graviton propagator. 



^.i/Apu 



{x; x') 



D-2 



iAcf(x; x' 



rpcf 

a/^^ pa 



iAci{x; x) . 
(129) 
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Table 6: Contractions from the zAcf part of the graviton propagator. 
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o 

I 


3 


b 
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3 


2 


a 


-4(IJr2) x) t[S] (x; x^T^ 


3 


2 


b 


x') dHA,i{x; x)} 


3 


2 


c 


-lK^^'{i[S]ix;x') ^tA,{{x;x)} 


3 


3 


a 


{%^)KH[S]{x;x')d-d'iA,i{x;x') 


3 


3 


b 


-U^)^'^^'[S]i^-^^W'^ciix;x) 


3 


3 


c 


+ ^K^Y^[S]{x;x')d'^^iA,,{x;x) 



We carry out the reduction in three stages. In the first stage the conformal 
part ll29l of the graviton propagator is substituted into the generic results from 
Table [5] and the contractions are performed. We also make use of gamma 
matrix identities such as Equation [32] and, 



S 



{x;x')^^ = {D-2)i\s]{x;x') and -f^J""" = --{D-l)-f^ 



^ (130) 

Finally, we employ relation [122l whenever c/ acts upon the fermion propagator. 
However, we do not at this stage act any other derivatives. The results of 
these reductions are summarized in Table El Because the conformal tensor 
factor [apTpl.] contains three distinct terms, and because the factors of 7° J^^^ 
in Table [5] can contribute different terms with a distinct structure, we have 
sometimes broken up the result for a given vertex pair into parts. These 
parts are distinguished in Tableland subsequently by subscripts taken from 
the lower case Latin letters. 
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In the second stage we substitute the fermion and conformal propagators, 



S 



"2vrf Ax« ' ^^^^^ 



T(^-l)(aa'V-T 



.A,(x;x')^^^^. (132) 

At this stage we take advantage of the curious consequence of the automatic 
subtraction of dimension regularization that any dimension-dependent power 
of zero is discarded, 

hm iAc{{x;x') = and hm d' iAc{{x; x') = . (133) 

In the final stage we act the derivatives. These can act upon the conformal 
coordinate separation Ax^ = x^ — x'^, or upon the factor of (aa')^~~ from 
the conformal propagator. We quote separate results for the cases where 
all derivatives act upon the conformal coordinate separation (Table [7]) and 
the case where one or more of the derivatives acts upon the scale factors 
(Table [8]). In the former case the final result must in each case take the form 
of a pure number times the universal factor. 



{aa'y °^^Ax^ 



(134) 



The sum of all terms in Table [7] is, 

-^i.,,, ^ ^r(§)r(f l)(-2Z,^+5Z,^)(Z,-l)(aa')'-*^ . 

(135) 

If one simply omits the factor of {aa'Y~~ the result is the same as in flat 
space. Although Equation 11351 is well defined for x'^^x^ we must remember 
that [S](a;; x') will be used inside an integral in the quantum-corrected Dirac 
equation shown by Equation [2H For that purpose the singularity at x'^ = x^ 
is cubicly divergent in D = 4 dimensions. To renormalize this divergence we 
extract derivatives with respect to the coordinate x^, which can of course be 
taken outside the integral in Equation [21] to give a less singular integrand. 



l^Ax^ \ l_ 

Ax2^ 2(D-l)\Ax2^' 



(136) 
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Table 7: Conformal iAct terms in which all derivatives act upon Ax^(x;x'). 



All contributions are multiplied by r(y)r(-j — l)(aa')^ ?. 
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(137) 
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Expression 11381 is integrable in four dimensions and we could take D = 4 
except for the explicit factor of 1/(D — 4). Of course that is how ultraviolet 
divergences manifest in dimensional regularization. We can segregate the 
divergence on a local term by employing a simple representation for a delta 
function, 



D-4 VAa;2^-6/ D-4[Ax2^-6 j r(f-l) D-4 



92 f ln(/i2Ax2) 1 z47rT/i^-4 ^^(x-x') , , 
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The final result for Table [7] is 

.2 1 



IK, 



+ 0(D-4) 



J ln(/i^Ax^) 

(2D2-5D+4)(aa')^~^ 
(D-2)(D-3)(D-4) 



(141) 



When one or more derivative acts upon the scale factors a bewildering 
variety of spacetime and gamma matrix structures result. For example, the 
[3— 2]b term gives, 
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(142) 
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The first term of Equation 11431 originates from both derivatives acting 
on the conformal coordinate separation. It belongs in Table The next 
three terms come from a single derivative acting on a scale factor, and the 
final term in Equation 11431 derives from both derivatives acting upon scale 
factors. These last four terms belong in Table [81 They can be expressed as 
dimensionless functions of D, a and a' times three basic terms, 
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■.(144) 



These three terms turn out to be all we need, although intermediate 
expressions sometimes show other kinds. An example is the [3—1] term. 
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Table S 
factors. 



Conformal iAd terms in which some derivatives act upon scale 
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All contributions are multiplied by 
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1, , a'i77° 



4 



'Ax2^-2 

Ax2^-2 



(146) 



As before, the first term in Equation 11461 belongs in Table [7l The second 
and third terms are of a type we encountered in Equation 11431 but the final 
term is not. However, it is simple to bring this term to standard form by 
anti-commuting the through either 7°, 



aa'H^-f°Y^x^j° = -aa'H^-f^'Ax^ - 2aa'H^Ar]j^ , 
= -aa'H^Y^x^ - 2(a-a')i^7° . 

Note our use of the identity {a — a') = aa'HAi]. 



(147) 
(148) 
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When all terms in Table [S] are summed it emerges that a factor of H^aa' 
can be extracted, 



[x-x') = —^TH-]iaa'y--{ (D^-7D-,.,,. , , 

+\{D-2)[a-a') x ^^-i(Z^-2)(3D-2)(a-a') x ^^^},(149) 



X 



-i(--)x^-i(^-)(3Z^-2)x^). ,150) 

Note the fact that this expression is odd under interchange of and x'^. 
Although individual contributions to the last two columns of Tableware not 
odd under interchange, their sum always produces a factor of a—a' = aa'HAr) 
which makes Equation 11501 odd. 

Expression 11501 can be simplified using the differential identities, 

Ar7^7^Ax^ _ / 7^Ax^ n 

Aa;2^ ~ 4(D-2)(D-1) VAx2^-4/ 

1 7^Ax^ , 1 ^^At] ^^^^^ 



2{D-1) Aa;2^-2 D~l Ax^^-^ ' 
j°Ar] -f^do f 1 



Ax2^-2 2{D-2) \Ax^^' 
The result is, 



(152) 



16(D-1) ' "VAx2^-4y 32V D-1 ^ "VAs2^-4 

We now exploit partial integration identities of the same type as those pre- 
viously used for Table [71 

^'^^^ = ( ^ ) = J-(J-) + 0(D-4) (154) 

Ax2^-4 2(Z}-3) VAx2^-6^ 2VAx2^ ^ ^' ^ ' 

7'^Ax^ -^d^ /I 



Ax2^-2 4(D-2)(D-3)(D-4) VAx2^-6 
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1 / 1 



(D-3)(D-4) ' 



(155) 



Ax2^-4 2(D-3)(L'-4) VAx2^-6 



(156) 



4 \ Ax2 / ' ""^"^ ^' ' r(f -1) (D-3)(D-4) ■ 
It is also useful to convert temporal derivatives to spatial ones using, 

^o^o =^ -? and dl = V' - 52 . (157) 
Substituting these relations in Equation 11531 gives, 

= — <^ -(D3-13D2 + 27D-12 

^ J 297r#(Z}-l)(Z}-3)(D-4) 1 ^ 

[6 ^92-2^92] ( 



— z 



-2D(D-2)^U^(x-x') + 



29-3-7r4 



'ln(;U^Ax^ 



Ax2 



+5^(V2-a2)(^)|+0(D-4) . (158) 



3.4 Sub-Leading Contributions from i5/\ 



A 



In this subsection we work out the contribution from substituting the residual 
A-type part of the graviton propagator in Table [5l 



D-3 



pa 



i8/\A{x]x') . (159) 



As with the conformal contributions of the previous section we first make 
the requisite contractions and then act the derivatives. The result of this 
first step is summarized in Table [91 We have sometimes broken the result 
for a single vertex pair into as many as five terms because the three different 
tensors in Equation ll59l can make distinct contributions, and because distinct 
contributions also come from breaking up factors of 7"J^'^. These distinct 
contributions are labeled by subscripts a, 6, c, etc. We have tried to arrange 
them so that terms closer to the beginning of the alphabet have fewer purely 
spatial derivatives. 
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Table 9: Contractions from the i6AA part of the graviton propagator 
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The next step is to act the derivatives and it is of course necessary to 
have an expression for i6AA{x;x') at this stage. From Equation [551 one can 
infer, 

2(5A^(x; x') = 

r(f + l) (aa')2-^ if^-2r(D-l)f /TT N , , ,1 
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Table 10: Residual i5/S.A terms giving both powers of Ax^. The 
two coefficients are Ai = ^^r(Y+i)x T{^){aa'y^~ and A2 = 



Function 


Vertex Pair 2-1 


Vertex Pair 2-2 
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(D-2)(D-3)2(D-4) 
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^D^2 ^ 1 1 r(n + D-l) 1 r(n+f -H) ^l/^n-f +2) 



(47r)f^Jn r(n+f) + 2 r(n + 2) V4 

In D = 4 the most singular contributions to Equation 11201 have the form, 
i5A^/Aa;^. Because the infinite series terms in Equation 11601 go like positive 
powers of Ax^ these terms make integrable contributions to the quantum- 
corrected Dirac equation in Equation [2H We can therefore take D = 4 for 
those terms, at which point all the infinite series terms drop. Hence it is only 
necessary to keep the ffist line of Equation 11601 and that is all we shall ever 
use. 

The contributions from iSAa are more complicated than those from zAcf 
for several reasons. The fact that there is a second series in Equation [160] oc- 
casions our Table [TOl These contributions are distinguished by all derivatives 
acting upon the conformal coordinate separation and by both series making 
nonzero contributions. Because these terms are special we shall explicitly 
carry out the reduction of the 2—2 contribution. All three 2—2 contractions on 
Table |9] can be expressed as a certain tensor contracted into a generic form, 

2 1 r 

SijSk£ + SikSje-jj—^Sie6jk\ x —daj\^i6AA{x;x')dki[S]{x;x')-fij . (161) 
So we may as well work out the generic term and then do the contractions 
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at the end. Substituting the fermion propagator brings this generic term to 
the form, 

Generic = —di'yj^i6AA{x; x')dki[S]{x; x')'yej , (162) 
''''^^%^j{i5A4x;x')d,{^^)^,} . (163) 



»7r 2 

Now recall that there are two sorts of terms in the only part of i6AA{x; x') 
that can make a nonzero contribution for D — A, 

^5AA,i.■y) . (164) 

z5Aa2{x;x') ^ ^Q^{-7rcot(fD)+ln(aa')}. (165) 

[Air) 2 i [-2) ' 

Because all the derivatives are spatial we can pass the scale factors outside 
to obtain, 
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= ^S^r(i?-2){-7rcot(|i^)+ln(aa')}9.7A^7.(^) ■ (168) 
To complete the reduction of the first generic term we note, 
1 ^ (t^\ ^ Ik Di^Ax^Ax^ 
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Hence the first generic term is, 
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Now we contract the tensor pref actor of Equation ll61l into the appropriate 
spinor-differential operators. For the first generic term this is, 
2 



+ 



a%7i^- (173) 



-VV+2V27°5o 



This term can be simphfied using the identities, 

liH = -7.7.^-2^= (1^-1)^-2^= (D-3)^ + 27% 

Applying these identities gives, 
2 



(174) 
(175) 
(176) 
(177) 



D-3 



2D 



D-A 
D-3 



D-3 

For the second generic term the relevant contraction is, 

2 



D-3 



d^P . (178) 



SijSke + SikSje- 



D-3 



T^ieSjk 



X dajdk p'ji 

D-5^^^^ 



^D-3^ 
D- 



VV-4(§^)V^ 



(179) 

D-3^ • ^ -^D-3^ • ^ • ^'''^ 

In summing the contributions from Table [10] it is best to take advantage 
of cancellations between Ai and A2 terms. These occur between the 2nd and 
3rd terms in the second column, the 4th and 5th terms of the 3rd column, 
and the 6th and 7th terms of the 3rd column. In each of these cases the 
result is finite; and it actually vanishes in the final case! Only the first term 
of column 2 and the 2nd term of column 3 contribute divergences. The result 
for the three contributions from [2—1] in Table [10] is, 

K^H'fi^-^ (Z}-i)r(f +1) 



2%f 



(D-3)2(D-4) 



aa 



x — x 



2 ^ 



ln(/i^Ax^ 
Ax2 



,-r4 + 21n(ii72^x2 



Ax2 



-0(D-4).(181) 
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Table 11: Residual iSAa terms in which all derivatives act upon Ax^(x;x'). 
All contributions are multiplied by ^-^^T{^+i)T{^){aa'Y~^ . 



I 


J 


sub 






1 A£'||27MAzp 


||Ax||27^Axi 








Az^-D 


2 


3 


a 




2D 








2 


3 


b 





1 


2 


-2d 


2 


3 


c 







D-3 


2D 
D-3 


3 


1 


a 


4(i:>-i)(i:>-2) 











3 


1 


b 












3 


2 


a 





4(ni) 

^ JJ — 6 ' 








3 


2 


b 




/D+lN 




-4(§^) 


3 


2 


c 




-|(D+l) 


-(D-l) 


2(d-i) 


3 


2 


d 








2(d-i) 


3 


3 


a 


^(D-l)(D-2) 
^ (D-3) 











3 


3 


b 












3 


3 


c 


\D-?,I 










3 


3 


d 


(D-l)(D-5) 




(Z3_5)(D_2) 


(D-5)(D-2) 


4(D-3) 




4(D-3) 


2(D-3) 


3 


3 


e 








-1(^-2) 



The result for the five contributions from [2—2] in Table [TO] is, 

^2^2^D-4 r(f + l) 



25^f (D-3)2(D-4) 



i\2-- 



+ 



2%4 





pln(/x2Aa;2)^ 











2 + ln(ii72/^a;2 



Ax2 



0(D-4). (182) 



As might be expected from the similarities in their reductions, these two 
terms combine together nicely in the total for Table [T0| 



[x: X 



r\2- 



2%f (D-SnD-A) 



■+2 



-{D-l)^+^]5^{x-x') 
ln(/i2Ax^)-| 



A,t2 
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Ax2 



0(D-4).(183) 



The next class is comprised of terms in which only the first series of iSA^ 
makes a nonzero contribution when all derivatives act upon the conformal 
coordinate separation. The results for this class of terms are summarized 
in Table [Til In reducing these terms the following derivatives occur many 
times, 



diiSAA^x; x) 
9oi5A^(x; x') 



, £r(- + lj(aa; 

:^r(§+i)(.af-« 



Ax' 



Ax^-2 

At] 

Ax^-2"2Ax^-4 

rD-2 



-d[i6AAix;x'), (184) 
aH 1 



9qz5A^(x; x') 



-^r(— + lj(aa 
Stt 2 ^ 



ij^-^r(D-i) 

^2^7rf r(f) 
At] a'H 
' Ax^-^~ 2Ax^-^ 
H^-^T{D-1) 



aH , (185) 



2^7rf 



r(f) 



a . (186) 



We also make use of a number of gamma matrix identities, 



7*7^Ax^7* 
7^AxS''Ax^7^Ax^' 



-D and 7V = -(^-1), 
{D-2)Y and 7^7* = (£'-1)7'- 
-As^ and (yAx*)^ = -||Aff , 
(D-l)7^Ax^-27*Ax* , 
||Af|| VAx/, - 2||Ax|| VAx^ • 



(187) 
2r , (188) 
(189) 
(190) 
(191) 



In summing the many terms of Table fTTl the constant K = D — 2/{D—3) 
occurs suspiciously often, 



^r(-+i)r(-)(aa')^ 



X • 



-2(D-1) 



D-1 



K 



7^ Ax 



V 2 



Ax2^-2 



+ 



iD-2) 



K 



YAxi 

Ax2^-2 



/D-2x ||Af||VAx^ {D-2){D-A) \\AxfYAxi 



Aa; 



2D 



P-3) 



Ax2^ 



..(192) 
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The last two terms can be reduced using the identities, 



|Af||27^A 



1 7^ Ax, 



Ax2^ 



+ 



yAxi 



7^ Ax 



2Ax2^-2 D-lAx2^-2 4(D-2)(D-1) VAx2^5-4 
D + lx 7^Axi , / 7* Ax 



if 



+ 



Ax2« 2VD-l^Ax2^-2 4(D-2)(D-1) VAx2^-4 
Substituting these in Equation 11921 gives, 



^ ^ (193) 
(194) 



x; X 



267r^ 



7^Axu 



Ax 



2D-2 



(D-2)(D2 



2 



1) + 



-10) 



2{D 
DK 1 yAx 



2(D-l)(D-3) 4(D-1) 
KV'^ ( 7'^Ax^ ^ , (L)-4)V2 



16(D- 



Ax2^- 



+ 



Ax2^-2 

7*Axj 



4(D-l)(D-3) Ax2^-4 



(195) 



We then apply the same formalism as in the previous sub-section to partially 
integrate, extract the local divergences and take D = 4 for the remaining, 
integrable and ultraviolet finite nonlocal terms. 



—I 



' DK 



x; X 



_D 

■ 2 



277rf 



-2(D-1) ^ 



P_3)P_4) 
M (L)-2)(D2- 



-5L> + 10) 



m2g2 
29-3-7r4 



-15^9^-4^(92] ( 



4(D-1) 2(D-l)(D-3) 
ln(/i2Ax2)x 1 



5^(x-x'l 



Ax2 



Ax^ 



-0(D-4). (196) 



The final class is comprised of terms in which one or more derivatives act 
upon a scale factor. Within this class we report contributions from the first 
series in Table [12] and contributions from the second series in Table [131 Each 
nonzero entry in the 4th and 5th columns of Table [12] diverges logarithmically 
like l/Ax^^"^. However, the sum in each case results in an additional factor 
of a — a' = aa'HAri which makes the contribution from Table [T2] integrable. 



D 
~2 



r(-+i)r 



3 / Ax2«-4 



1 /3D-4n 7^Axi7'^Ax^7°Ar/ (L)-l)(D-4) 7^Ax^ 
2V D-Z ' Ax2^-2 + 4(D-3) Ax2^-4 



(197) 
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Table 12: Residual iSAa terms in which some derivatives act upon the scale 
factors of the first series. The factor ^^^^ V{^+i)V{^){aa'Y~~ multiplies all 
contributions. 



T 
i 


T 
J 


sub 


//70 










Ax'20-4 


2 


1 




2(§e|)a' 


my 





2 


3 


a 


-my 


my 





3 


1 


a 








(D-l)(D-4) 
2{D-3) 


3 


1 


b 





my 





3 


2 


a 


-my 


-2(fei)a 





3 


3 


a 








(D-l)(D-4) 
4(D-3) 


3 


3 


b 










3 


3 


c 











Table 13: Residual i5^A terms in which some derivatives act upon the 
scale factors of the second series. All contributions are multiplied by 



I 


J 


sub 




AxO+2 


Q2/ H7O N 

" IaxO-2'' 


2 


1 




-^my 


-my 





2 


3 


a 


my 


my 





3 


1 


a 








(D-l)a 
(D-2)(D-3) 


3 


2 


a 


m)- 








This is another example of the fact that the self-energy is odd under inter- 
change of and x'^. 

The same thing happens with the contribution from Table | 



[ tic > tic 



m - Daa'l (-^) 



D \ 7* Axj7^Ax^7°Ar7 



D-3^ Ax^ 



Ax^^ 



Ha 



>.(198) 
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We can therefore set D = A, at which point the two Tables cancel except for 
the delta function term, 



(199) 

It is worth commenting that this term violates the reflection symmetry of 
Equation 11031 In D = 4 it cancels the similar term in Equation 11271 



3.5 Sub-Leading Contributions from i6A 



B 



In this subsection we work out the contribution from substituting the residual 
5-type part of the graviton propagator in Table 0, 



Q/jApcr 



K^'^pp + OTlp^ + + zSAb . (200) 



As in the two previous sub-sections we first make the requisite contractions 
and then act the derivatives. The result of this first step is summarized in 
Table [HI We have sometimes broken the result for a single vertex pair into 
parts because the four different tensors in (12001) can make distinct contribu- 
tions, and because distinct contributions also come from breaking up factors 
of 7"J^^. These distinct contributions are labeled by subscripts a, b, c, etc. 

2(5Ab(x; x') is the residual of the S-type propagator of Equation [8^ after 
the conformal contribution has been subtracted, 



iSAB(x; x) 



H^T{^) {aa'y-'^ g^-^ r(D-2) 
IGttt (4^)f r(f) 

^ r r(n+f ) .yx „-f +2 T{n + D~2) ^y^ 



■ (4vr)f r(n + 2) r(n+f) j" ^'"'^ 

As was the case for the i6AA{x] x') contributions considered in the previ- 
ous sub-section, this diagram is not sufficiently singular for the infinite series 
terms from i6AB{x\ x') to make a nonzero contribution in the D = A limit. 
Unlike i5AA{x\ x'), even the n = terms of iSAsix; x') vanish for D = 4. This 
means they can only contribute when multiplied by a divergence. 

Contributions from the [2—2] vertex pair require special treatment to take 
advantage of the cancelation between the two series. We will work out the 
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Table 14: Contractions from the idAs part of the graviton propagator. 



I 


J 


sub 


■T rCX.3 / \ • r Ol / /\ tO<J f l\ r rrifil 'fA / /\ 

'^(a;) x') iVj (x') [a/jT,^] «5Aij(x; x') 


2 


1 







2 


2 


a 


-|«:2a;{7Wz[5](x;x')7fc^5AB(a;;a;')} 


2 


2 


b 




2 


3 


a 


-lK,^'ykdoi[S]{x; x')'y'' d'^iSAB^x; x') 


2 


3 


b 


\K'^-1^d'^i5AB{x; x') dki[S]{x] x')-i^ 


2 


3 


c 


-^K'^^y'^dkiSAsix; x') dQi[S]{x] x')7° 


2 


3 


d 


1^2709^-2 [^] {x- a;')7° dki5AB{x; x') 


3 


1 







3 


2 


a 


|K29^{7'=i[S'] (x; x')7/fc 9o25Ab(x; x')} 


3 


2 


b 


i/^27'=9U^[5](x; x')7° d,i5AB{x- x')} 


3 


2 


c 


-\K'^-i^d'Q{i[S]{x; x')-i^ dkiSABix] x')} 


3 


2 


d 


-lK^dk{i^t[S]{x; x')7° dH6AB{x; x')} 


3 


3 


a 


-j^K^'yki[S] (x; x')7''^o5o^'5Aij(x; x') 


3 


3 


b 


3^«:V^[5](x; x')7' dkd'oi6AB{x; x') 


3 


3 


c 


-^K^-f''i[S]{x; x')7° dodkiSAsix; x') 


3 


3 


d 


^K^-f^i[S]{x; x')7°V2z5Ab(x; x') 



a" term from Table [T] 

^'^'r(f-i)^, 



2-2 



IGtt 2 



9^ z5Ab(x;x') 7V^-f ^f^c 



1 



Ax^-2 



^'^'r(f-i)^, 



IGtt 2 



'^d'Ji6AB{x;x')[-3do^+-f'V^ + {D-l)^d, 



Ax^' 



(202) 
■.(203) 



A key identity for reducing the [2 — 2] terms involves commuting two deriva- 
tives through l/Ax-^~^, 



(204) 



Ax^-4 " 

This can be used to extract the derivatives from the first term of i6AB{x; x'), 
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Table 15: Residual iSAs terms in which all derivatives act upon Ax'^{x]x'). 



All contributions are multiplied by ^^r2(f )(£)-4)(aa')^^ t. 



1 


J 


sub 






IIAJII--'' A,r„ 

1 1 — 1 1 / — /' 












2 


3 


a 


(D-l)2 


-{D+1) 


-D(D-l) 


2D 


2 


3 


b 


(D-l) 


-2D+1 


-D 


2D 


2 


3 


c 





-{D-l) 


-D 


2D 


2 


3 


d 





-1 


-D 


2D 


3 


2 


a 


-2(D-l)(D-2) 


3D-5 


2(D-1)2 


-4(D-1) 


3 


2 


b 


-(D-l) 


3(D-1) 


2(D-1) 


-4(D-1) 


3 


2 


c 





2D-3 


2(D-1) 


-4(D-1) 


3 


2 


(1 




2i)-i 


2U)-1) 


-J-U-'-i) 


3 


3 


a 


i(D-l)(D-3) 


-(D-3) 


-i(D-l)(D-2) 


(D-2) 


3 


3 


b 







-5 (-0-2) 


(D-2) 


3 


3 


c 







-|(i?-2) 


(D-2) 


3 


3 


d 


i(D-l) 


-(D-l) 


-|p-2) 


(D-2) 


Total 






-i(D-l)(D-2) 


3(D-2) 


i(D+2)(D-2) 


-4(D-2) 



at which point the result is integrable and we can take = 4, 
2-2 



2%^ 



(aa'Y-T 



Ax 



D-4 



Ax^- 



m^g^r(f+i)r(f-i) 

2%^ D-3 



(aa') 



/\2- 



X 



- - {D-4)Ha') (-30*0 ^+^'V'+{D - 1) ^do 
■7%(39o^ + V^)[-^] +0(1^-4). 



Ax^^- 



2%4 



.,(205) 



(206) 
(207) 



Of course the second term of ISAb is constant so the derivatives are already 
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extracted, 



2-2 



2%4 



-7 



'5o(392 + V 



1 



Ax2 



0(D-4) . 



,(208) 
(209) 



Hence the total for [2— 2]^ is zero in D = 4 dimensions! 

The analogous result for the initial reduction of the other [2—2] term is, 



2-2 



m2r(f-i) 

167r~ 



X 



(210) 



The results for each of the two terms of iSAb are, 



2-2 



1 m2ff2r(f +i)r(f-i 



2%^ 



-27%V2+^V'+^ao') 



1 



Ax2^-6 



(-27% V' +^V' +^5o' 



2-2 



2^4 
iifH^lXD-2 
2^+%^ D-2 



Aa;2 

27°9oV2-?V2 " 



+ 0(D-4) 
1 



Ax^-2 



2%^ 



(27W -^9o^) [^] + OiD-4) . 



(211) 
(212) 
(213) 
(214) 



Hence the entire contribution from [2—2] vanishes in D = A. 

The lower vertex pairs all involve at least one derivative of iSAs, 



diiSAsix; x') 
doiSAB^x; x ) 
&Qi5AB{x] x') 



IGvr 2 

^!n#*(fl-4)(aa')-* 

lOTT 2 

--^(D-4)(aa')^-- 
lovr 2 



=-5:5A5(x;x'),(215) 

[ Ar] aH \ 

jAx^-2~2Ax^-4J' 

f A77 a'i/ 1 



Ax^-2 2Ax^-4 r- 



(216) 
(217) 



These reductions are very similar to those of the analogous iSA^ terms. We 
make use of the same gamma matrix identities of Equations 1187111911 that 
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were used in the previous sub-section. The only really new feature is that 
one sometimes encounters factors of /S.7f which we always resolve as, 



-Ax^ + IIAf I 



(218) 



Table [15] gives our results for the most singular contributions, those in which 
all derivatives act upon the conformal coordinate separation Ax^. 

The only really unexpected thing about Table [15] is the overall factor of 
(D — 2) common to each of the four sums. 



X] X 



t\-){D-2){D-A){ 



2% 



D 



+3 



2 
1 



aa 



J\2- 



D_ 



1 'y'^Ax 
-(D-1) ^ " 
2^ ^Ax2^-2 



'j^Axi 



I Axp7'^Aa;^ ||Axp7*Axj 



Ax2^ 



Ax2^ 



.(219) 



As with the result of Table [TT] we use the differential identities 1193111941 to 
prepare the last two terms for partial integration. 



[^){D-2){D-4){aar-- 




16 (D-l)(D-3)'^ 2{D-l){D-3) 
The expression is now integrable so we can take D = 4 



j Ax2^^ 



6 ■ 



(220) 



— t 



[x: X 



2%^ 



(221) 



Unlike the iSA^ terms there is no net contribution when one or more 
of the derivatives acts upon a scale factor. If both derivatives act on scale 
factors the result is integrable in D = 4 dimensions, and vanishes owing to 
the factor of (_D — 4)^ from differentiating both a^~~ and a'^~~. If a single 
derivative acts upon a scale factor, the result is a factor of either {D — 4)a 
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or (D— 4)a' times a term which is logarithmically divergent and even under 
interchange of and x'^. As we have by now seen many times, the sum of 
all such terms contrives to obey reflection symmetry of Equation 11031 by the 
separate extra factors of (D — 4)a and {D — 4)a' combining to give, 

{D-A){a - a!) = {D-A)aaHAr] . (222) 

Of course this makes the sum integrable in D = 4 dimensions, at which point 
we can take D = 4 and the result vanishes on account of the overall factor of 
(D-A). 



3.6 Sub-Leading Contributions from i6A 



c 



The point of this subsection is to compute the contribution from replacing 
the graviton propagator in Table [5] by its residual C-type part. 



C- 



{D-2){D-3) + + ^D^) 

(223) 

As in the previous sub-sections we first make the requisite contractions and 
then act the derivatives. The result of this first step is summarized in Ta- 
ble [161 We have sometimes broken the result for a single vertex pair into 
parts because the four different tensors in Equation 12231 can make distinct 
contributions, and because distinct contributions also come from breaking 
up factors of 7"J^'^. These distinct contributions are labeled by subscripts 
a, b, c, etc. 

Here i6Ac{x] x') is the residual of the C-type propagator of Equation 
after the conformal contribution has been subtracted, 



t5Ac(x;x') = rrf 3)r( 1) \ ' ^ + ^ 



T{n + D-3) 




(1)1.(224, 



,4/ J 

As with the contributions from i5AB{x\x') considered in the previous sub- 
section, the only way i6Ac{x\x') can give a nonzero contribution in D = 4 
dimensions is for it to multiply a singular term. That means only the n = 
term can possibly contribute. Even for the n = term, both derivatives must 
act upon a Aa;^ to make a nonzero contribution in D = 4 dimensions. 
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Table 16: Contractions from the i5/S.c part of the graviton propagator. 



I 


J 


sub 


iVT^ix) %\S\{x\ x'\ %V^^''{x'^ LflT^l i^Arfx; x'\ 


2 


1 

X 


o 


(D-3)(D-2) ^ X ) lUl^c , X J 


2 


1 


b 


J n^d' {i^dni\S]{x\x')Y i5Ar{x\x')\ 


2 


2 


a 


d 5^{x—x') iSAcix; x) 

2(D— 3)(D— 2) " ^ ' ov ) / 


2 


2 


b 


~2(n si^^7°^o ^^(x-a;') i(5Ac(a;; x) 


2 


2 


c 


+ 9rn ,,/t'<9'{7%i[5](x;x')7^^5Ac(a:;a:')} 

2( U — o) fi ' ^ I i \ 1 / I ^ \ 1 J J 


2 


2 


d 


^ ( n"3 ) /^'^o {7°^o i [-5] (x; a;')7° i(5Ac (x; x') } 

2 ^ U — O UL./ L JV ' // "^V ' /J 


2 


3 


a 


.rr.^'^T}! ..K^5^{x x') d' iSArix: x') 

4(D—3)(D—2) ^ I T »o V"*^ ) / 


2 


3 


b 


+ 4FTflT'^' ^''{x-x')-i%if>Ac{x;x') 


2 


3 


c 


+ \{n'\)i^^l'^doi\S]{x; x') d'iSAcix; x') 


2 


3 


d 


lin~l)'^^l^doi\S](x; x'Wd'i6Ac(x; x') 


3 


1 


a 


^ K^a' WiSAr(x: x') i\S](x: x'h^'} 


3 


1 


b 


+ 2in^^%b% ^^^c{x; x') i[S] {x; x')^} 


3 


2 


a 


4iD%XD-2)^"%{^i^^c{x-, x') i[S\ (x; x')^} 


3 


2 


b 


-\{%^W&,{^i5Ac{x- x') i[S]{x- x')7°} 


3 


2 


c 


i^^cix; x') i[S] (x- x')r} 


3 


2 


d 


mM)^'9',{Yd, tSAcix; x') t[S]ix; x'h^} 


3 


3 


a 


^^S£g^K'Y^[S]{x■,x')^,i^'^6Ac{x■,x') 


3 


3 


b 


-|(§E|)'^'7'^^[^](a:; x')d,^W'j i5Ac{x; x') 


3 


3 


c 


--siM)'^'y^[S]ix; x')d, ^ i5Ac{x- x') 


3 


3 


d 


+|(§^)'^V«['5] {x- x')da'd'^ i5Ac{x- x') 
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Those of the [2—1] and [2—2] vertex pairs which are not proportional to 
delta functions after the initial contraction of Table [12] all contrive to give 
delta functions in the end. This happens through the same key identity 12041 
which was used to reduce the analogous terms in the previous subsection. In 
each case we have finite constants times different contractions of the following 
tensor function, 



d'Ji6Ac{x;x')dad[3 



r(D-3) 

{An) 2 i [^) 



+ 



IGtt 2 ^ 2 



V 2 



Ax 



D-A 



1 



Ax 



D-2 



(225) 



H^~' T{D-3) ^ 
(47r) 



Ax 



D-2 



x(a;-i(B-4)ffa' 



4(D-3) 4{D-3) 



(226) 



167r2 ^ 



Ax"' 



j^2 



1 



Ax2 



16 



Vapd^S\x-x') + 0{D-A) . 



Ax2-D-6J 

0{D-A), (227) 
(228) 



It remains to multiply Equation 



take the appropriate contraction. For example, the [2—1];, contribution is 

-2 ,rtD 1^ ^^2 



K 



D-3 



A R 

Air 2 



by the appropriate prefactors and 
i, the [2— l]fe cor 

Vapd^S'^ix-x') 



16 



167r2 



X 



{x-x') +0{D-A) . (229) 



We have summarized the results in Table [T71 along with all terms for which 
the initial contractions of Table [16] produced delta functions. The sum of all 
such terms is, 



^'^\-l^-\¥y{x-x') + 0{D-4) . 



167r2 L 8 ' 4 

All the lower vertex pairs involve one or more derivatives of iSAc, 



(230) 



ddSA 



c 



^2 '\D-6){D-4){aa'f' 



Ax' 



327rT 



Ax^-2 



-d[t5Ac , (231) 
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Table 17: Delta functions from the i5^c pai't of the graviton propagator. 
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sub 




^^5\x-x') 


2 


1 


a 


1 

2 





2 


1 


b 


1 
1 

4 





2 


2 


a 


1 

4 





2 


2 


b 


1 
2 


1 

2 


2 


2 


c 


1 
8 





2 


2 


d 


1 

4 


1 

4 


2 


3 


a 








2 


3 


b 








Total 






3 
8 


1 

4 



Table 18: Residual i^A^ terms in which all derivatives act upon Ax^(a;; x'). 
All contributions are multiplied by |g^r(f )r(f -l)^^=gff=^(aa')^~^. 
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sub 






||a;|p7'^Ax^ 


||a;|p7'Azi 




Ax2D-2 


Aa;2D 


Aa;2D 


2 


3 


c 


-(D-l)2 


D(D-l) 








2 


3 


d 





(Z)-l)(Z?-2) 


D(D-2) 


-2D{D-2) 


3 


1 


a 


4(D-1) 











3 


1 


b 


-2(D-1) 


-2(D-4) 








3 


2 


a 


-2(D-1) 











3 


2 


b 


2(D-l)(D-2) 


-2(D-l)(D-2) 








3 


2 


c 




(D-4) 








3 


2 


d 





-(2D-3)(D-2) 


-2(D-l)(D-2) 


4(D-l)(D-2) 


3 


3 


a 













3 


3 


1) 




i(l)-l)U)-i) 








3 


3 


c 




i(D-l)(D-4) 








3 


3 


d 


-i(D-l)(D-2) 


(D-2) 


i(D-2)2 


-(D-2)2 


Total 






i(D-l)(D-2) 


2(i:)-l)-D(D-2) 


-i(D-2)2 


(D-2)2 
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■ D 



327r- 



, (D-6)(D-4)(aa')--/^^^ ^ 



■ D 



327rT 



-{D-6){D-A){acL 



2AxD-\ 
At] a'H 



(232) 
■. (233) 



Their reduction follows the same pattern as in the previous two sub-sections. 
Table [T8] summarizes the results for the case in which all derivatives act upon 
the conformal coordinate separation Ax^. 

When summed, three of the columns of Table [TS] reveal a factor of (-D— 2) 
which we extract, 



—I 



(x: X 



(D-2)(D-4)(D-6) 



(D-3) 



D 



YAxi 



Ax2^-2 



1,^ JIAfllVAx. 



iD-2) 



\Ax\\^Y^Xi 
A^ 



(234) 



We partially integrate Equation 12341 with the aid of Equations 1193111941 and 
then take D = 4, just as we did for the sum of Table [T^ 



ix: X 



( D 7''Ax^ 



IK 



H'^.D^.D MD-2){D-A)iD-6) 



X 



4 Aa;2^-2 



28vr^ V2/ V2 
+ " 



{D-3) 



D-2^ 2{D-1) 



7*Axj 
Ax2^- 



^r(^)r(^-i 

~)3 



+ 



7*Axj 

(Z}-1) VAx2^-4V ' 4(Z}-1) VAx2^-4 

160^^' 



(D-l)(D-3)2 
D-4 



8(D-2) 



16 



Ax 



2D- 



2V {^^^^ + >1a^ + ^(^-4)- 



(235) 



(236) 
(237) 



As already explained, terms for which one or more derivative acts upon a 
scale factor make no contribution in D = 4 dimensions, so this is the final 
nonzero contribution. 
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3.7 Renormalized Result 

The regulated result we have worked so hard to compute derives from sum- 
ming expressions [EZl [Mil [ISH [Ha [M IMl [230] and [2371 



167r2 



X -31n(a)^5^(a;-x') 



2%4 



Ax2 



15 



ln(/i2Ax2 



In(ii72^x2) 



Ax2 



Ax2 



Ax2 

0(D-4). (238) 



The various D-dependent constants in Equation 12381 are. 

■^-4 r(f-i) r _ . 2 ^ 



&2 



287r# 


(D-3)(D- 


4) 


/.^-^ 


r(f+i) 




2%f 


(D-3)(D- 


4) 




r(f+i) 





-2D + 1- 



D-2 



24 6 



35 



D D-1 D-3 



-D + 3 + 



9 



2%f (D-3)(D-4) [ ' ~' ' D-3 
H^-^r{D-l)( (D + l)(D-l)(D-4) TT 
(4vr)f r(f) 



(239) 
(240) 
(241) 



/7tD 

X — cot 

2 V 2 



2(D-3) 

(D-1)(D3-8D2 + 23D-32) 7 



8{D-2y{D-3y 



- (4.)f r(f) iil^-DT 



TT /TlD 

, X — cot 

3^ 2 V 2 



+ 



3 {D^-6D- 



4(D-2)2(D-3)2 2 



(243) 



In obtaining these expressions we have always chosen to convert finite, D = 4 
terms with acting on l/Ax^, into delta functions, 



-Ax2- 



(244) 
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Figure 3: Contribution from counterterms. 



All such terms have then been included in 62 and 63. 

The local divergences in this expression are canceled by the BPHZ coun- 
terterms enumerated at the end of section 3. The generic diagram topology 
is depicted in Figure [3], and the analytic form is, 



—I 



-ictm 



iCiij 5 



D, 



X — X 



(245) 



1=1 



-K 



''[ai{aa)-^^d'^ + a2D{D-l)H^^ + a3H'^^}6^{x-x'). (246) 

In comparing Equation 12381 and Equation 12461 it would seem that the simplest 
choice for the coefficients a,- is, 



tti = Pi 



a2 



/32 + b2 



and eta = /^s+fts . 



(247) 



DiD-1] 

This choice absorbs all local constants but one is of course left with time 
dependent terms proportional to In(aa'), 

1 In(aa') 



Pi(aa'Y 2 — ai[aa 
p2{aa'f-^ +b2- D{D-l)a2 



aa' 

7.5 



0{D-A) , (248) 



267r2 
1 



\n{aa) +0{D-A) , (249) 



In(aa') + 0(D-4) . (250) 



Our final result for the renormalized self-energy is 

fln(aa') 15 



—I 



IK 



2 



2%4 



, aa 



2^%"^ \ aa' 



^ + — \n{aa')H^^-7\n{aa')H^^Y^{x~x') 

2 A^2\ 



ln(/i2Ax2)^ ^ ^2^2 



Ax2 

52 A av72 



15 



+ 8^a"-4^V^ 



4 [V 2 
ln(ii72^a;2) 



Ax2 



-7^V' 



ln{jj,^Ax 



Ax2 



Ax^ 



(251) 
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4 CORRECTING THE MODES 



It is worth summarizing the conventions used in computing the fermion self- 
energy. We worked on de Sitter background in conformal coordinates, 



a^{ri)(^—drf' + dx-dxj where a{ri) 



Hrj 



(252) 



We used dimensional regularization and obtained the self-energy for the con- 
formally re-scaled fermion field, 



(253) 



The local Lorentz gauge was fixed to allow an algebraic expression for the 
vierbein in terms of the metric [40| . The general coordinate gauge was fixed 
to make the tensor structure of the graviton propagator decouple from its 
spacetime dependence [H], [50] . The result we obtained is. 



(x; xt) ■ 

.2 



K 



2%^ 



[aa 



ln(/i2Ax2)^ ^ K^H^ 



Ax^ 



2H^ 



ln{\H^Ax^ 



Ax2 



'15 



ln(/i^Ax^ 



Ax2 



- Ax2 
+ 0(/t^).(254) 



where = IGirG is the loop counting parameter of quantum gravity. The 
various differential and spinor-differential operators are. 



(255) 



where rj'^'^ is the Lorentz metric and 7^^ are the gamma matrices. The con- 
formal coordinate interval is basically Ax^ = {x — x')^{x — x'Yri^j,, up to a 
subtlety about the imaginary part which will be explained shortly. 
The linearized, effective Dirac equation we will solve is. 



i^ij^j{x) - / rfV 



[x] x') '^j{x') 







(256) 



In judging the validity of this exercise it is important to answer five questions: 



1. How do solutions to Equation 12561 depend upon the finite parts of 
counterterms? 
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2. What is the imaginary part of Ax^? 



3. What can we do without the higher loop contributions to the fermion 
self-energy? 

4. What is the relation between the (D-number, effective field Equation 
12561 and the Heisenberg operator equations of Dirac + Einstein? and 

5. How do solutions to Equation 12561 change when different gauges are 
used? 

In next section we will comment on issues 1-3. Issues 4 and 5 are closely 
related, and require a lengthy digression that we have consigned to section 2 
of this chapter. 



4.1 The Linearized Effective Dirac Equation 

Dirac + Einstein is not perturbatively renormalizable [18J, so we could only 
obtain a finite result by absorbing divergences in the BPHZ sense [191 1201 
[2T1 [22] using three counterterms involving either higher derivatives or the 
curvature R = 12 H^, 

- ^^H^{jH^^9^ + a2D{D-l)^+a3'^}6''{x-x') . (257) 

No physical principle seems to fix the finite parts of these counterterms so 
any result which derives from their values is arbitrary. We chose to null local 
terms at the beginning of infiation (a = 1), but any other choice could have 
been made and would have affected the solution to Equation 12561 Hence 
there is no point in solving the equation exactly. However, each of the three 
counterterms is related to a term in Equation 12511 which carries a factor of 
In(aa'), 

In(aa') 



H^aa' ^ H'^ 



aa 



a2/^(£'-l)^^ y In(aa')^, (259) 

ag'p^ -71n(aa')?. (260) 

Unlike the a^'s, the numerical coefficients of the right hand terms are uniquely 
fixed and completely independent of renormalization. The factors of In(aa') 
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on these right hand terms mean that they dominate over any finite change in 
the Oj's at late times. It is in this late time regime that we can make reliable 
predictions about the effect of quantum gravitational corrections. 

The analysis we have just made is a standard feature of low energy ef- 
fective field theory, and has many distinguished antecedents [23l IMl l25l [261 
ETIEaESlEniEIlEaESlIMlESlESlEZI. Loops of massless particles make 
finite, nonanalytic contributions which cannot be changed by counterterms 
and which dominate the far infrared. Further, these effects must occur as 
well, with precisely the same numerical values, in whatever fundamental the- 
ory ultimately resolves the ultraviolet problems of quantum gravity. 

We must also clarify what is meant by the conformal coordinate interval 
Ax^(x; x') which appears in Equation l251[ The in-out effective field equations 
correspond to the replacement, 

Ax^(x; x') — > Ax^_|_(x; x') = \\x — x \\'^ — {\ r] — t]' \ —iS^ . (261) 

These equations govern the evolution of quantum fields under the assumption 
that the universe begins in free vacuum at asymptotically early times and 
ends up the same way at asymptotically late times. This is valid for scattering 
in fiat space but not for cosmological settings in which particle production 
prevents the in vacuum from evolving to the out vacuum. Persisting with 
the in-out effective field equations would result in quantum correction terms 
which are dominated by events from the infinite future! This is the correct 
answer to the question being asked, which is, "what must the field be in order 
to make the universe to evolve from in vacuum to out vacuum?" However, 
that question is not very relevant to any observation we can make. 

A more realistic question is, "what happens when the universe is released 
from a prepared state at some finite time and allowed to evolve as it will?" 
This sort of question can be answered using the Schwinger-Keldysh formalism 
[HIISIESIIZZIESIIIHIIMIEI]. Here we digress to briefiy derive it. To 
sketch the derivation, consider a real scalar field, (p{x) whose Lagrangian (not 
Lagrangian density) at time t is L[(f{t)]. The well-known functional integral 
expression for the matrix element of an operator Oi [ip] between states whose 
wave functionals are given at a starting time s and a last time i is 



T*(Oi\ 



The T*-ordering symbol in the matrix element indicates that the operator 
Oi [ip] is time-ordered, except that any derivatives are taken outside the time- 



64 



ordering. We can use Equation 12621 to obtain a similar expression for the 
matrix element of the antz-time-ordered product of some operator 02[v^] in 
the presence of the reversed states, 



T*(0l[ip])\<^)\ (263) 
[dp] 02[p] ^p{()] e-*/>*^['^Wl ■^*[ip{s)] . (264) 



Now note that summing over a complete set of states $ gives a delta 
functional, 



5 



^_(£)-^+(£) 



(265) 



Taking the product of Equation 12621 and Equation I264[ and using Equation 
12651 we obtain a functional integral expression for the expectation value of 
any anti-time-ordered operator O2 multiplied by any time-ordered operator 



r(o2[p])T*(oM)\^) 



[dp+][dpJ,6 p^{l)-p+{i) 



:b_]0,b+]^*b4^)]e^^»^'*{''""'^'"''^-'^^}^b4^)].(266) 



This is the fundamental relation between the canonical operator formalism 
and the functional integral formalism in the Schwinger-Keldysh formalism. 
The Feynman rules follow from Equation 12661 in close analogy to those for 



in-out matrix elements. Because the same field is represented by two different 
dummy functional variables, ip±{x), the endpoints of lines carry a ± polarity. 
External lines associated with the operator 02[</'] have — polarity whereas 
those associated with the operator Oi [ip] have + polarity. Interaction vertices 
are either all + or all — . Vertices with + polarity are the same as in the 
usual Feynman rules whereas vertices with the — polarity have an additional 

minus sign. Propagators can be ++, — h, H — and . 

The four propagators can be read off from the fundamental relation 12661 
when the free Lagrangian is substituted for the full one. It is useful to denote 
canonical expectation values in the free theory with a subscript 0. With this 
convention we see that the ++ propagator is just the ordinary Feynman 
propagator. 



T(p{x)p{x'] 



iA(x: x') 



(267) 
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The other cases are simple to read off and to relate to the Feynman propa- 
gator, 



^A__^_{x■, x') 
iA_,._(a;; x') 
iA__(x;x') 



= (fl 



(p{x)(p{x'] 



e{t~t')iA{x; x')+e{t'-t) \iA{x; x') 



ip{x')ip{x) 



e{t-t') \iA{x; x')] +e{t' -t)iA{x; x'), 



T[(fix)(p{x' 



iA{x; x') 



(268) 

(269) 
(270) 



Therefore we can get the four propagators of the Schwinger-Keldysh formal- 
ism from the Feynman propagator once that is known. 

Because external lines can be either -|- or — every N-point IPI function of 
the in-out formalism gives rise to 2^ IPI functions in the Schwinger-Keldysh 
formalism. For example, the IPl 2-point function of the in-out formalism — 
which is known as the self- mass-squared M'^{x; x') for our scalar example — 
generalizes to four self-niass-squared functions. 



M\x;x') 



Ml^{x-x'). 



(271) 



The first subscript denotes the polarity of the first position and the second 
subscript gives the polarity of the second position x'^. 

Recall that the in-out effective action is the generating functional of IPI 
functions. Hence its expansion in powers of the background field 4>{x) takes 
the form, 

T[(j)] = S[(t)] -\jd^x J d^x'(t){x)M\x; x')4){x') + 0(0^) , (272) 

where S\(p\ is the classical action. In contrast, the Schwinger-Keldysh effec- 
tive action must depend upon two fields — call them (p+{x) and (p-{x) — 
in order to access the different polarities. At lowest order in the weak field 
expansion we have, 

d^x [ d*x' I 0+ {x) MlJ,x- x')(t)+{x') + (t)+{x) Ml_{x; x')(p_ {x') 



+(/)_(x)iVr^(x;x')0+(a;') + 0_(a;)iVr_(a;;a;>_(a;')^+O(0^). (273) 
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d^x'M\x;x')(f){x') + 0{(f)' 



(274) 



The effective field equations of the in-out formahsm are obtained by vary- 
ing the in-out effective action, 

sm _ 5s[(i>] 

6(f){x) 6(j){x) 

Note that these equations are not causal in the sense that the integral over 
x'^ receives contributions from points to the future of x^. No initial value 
formalism is possible for these equations. Note also that even a Hermitian 
field operator such as (p{x) will not generally admit purely real effective field 
solutions (f){x) because IPI functions have imaginary parts. This makes the 
in-out effective field equations quite unsuitable for applications in cosmology. 
The Schwinger-Keldysh effective field equations are obtained by varying 



with respect to 



S(j)+{x) 



^ and then setting both fields equal, 

6S[<j>] 



4>±-- 



6(j){x) 



d^x' 



[X 



i+0( 



(275) 



The sum of M^^(a;; x') and M^_(x; x') is zero unless x'^ lies on or within the 
past light-cone of x^. So the Schwinger-Keldysh field equations admit a well- 
defined initial value formalism in spite of the fact that they are nonlocal. 
Note also that the sum of M^^(x; x') and M^_(x; x') is real, which neither 
IPI function is separately. 

From the preceding discussion we can infer these simple rules: 

• The linearized effective Dirac equation of the Schwinger-Keldysh for- 
malism takes the form Equation 12561 with the replacement. 



[x] X ) + jSj 



(276) 



The ++ fermion self-energy is Equation 12511 with the replacement Equa- 
tion [2611 and 

The +- fermion self-energy is. 



(277) 




with the replacement. 



Ax' 



X 



^l|2 

X 



- [r] -7]' + i5y 



(278) 
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The difference of the ++ and +- terms leads to zero contribution in Equation 
I unless the point x''^ lies on or within the past light-cone of x'^. 



We can only solve for the one loop corrections to the field because we lack 
the higher loop contributions to the self-energy. The general perturbative 
expansion takes the form, 



and 



£=1 



2e 



(x; x') 



(279) 



One substitutes these expansions into the effective Dirac equation in Equa- 
tion [256] and then segregates powers of k^, 

i^^\x) = , i^-^\x) = J d^x'[j:^]{x;x')^\x') et cetera. 

(280) 

We shall work out the late time limit of the one loop correction "^Kt], x; k, s) 
for a spatial plane wave of helicity s, 



^^i{r],x;k,s) 



-ikr) 



Ui{k, s)e 



ik-x 



where k^'^ljUj{k, s) = k'^^jUj{k, s) . 

(281) 



4.2 Effective Field Equations 

The purpose of this section is to elucidate the relation between the Heisenberg 
operators of Dirac + Einstein — i'iix), ipii^x) and h^^{x) — and the (D- 
number plane wave mode solutions \E'i(x; k, s) of the linearized, effective Dirac 
equation in Equation 12561 After explaining the relation we work out an 
example, at one loop order, in a simple scalar analogue model. Finally, we 
return to Dirac -|- Einstein to explain how "^{{x; k, s) changes with variations 
of the gauge. 

One solves the gauge-fixed Heisenberg operator equations perturbatively, 

h^Ax) = h1^,ix) + Khl,{x) + K^hl^{x) + ... , (282) 
^,{x) = ilj^{x)+K^l{x)+K^ilj^{x) + ... . (283) 

Because our state is released in free vacuum at t = [r] = —1/H), it makes 
sense to express the operator as a functional of the creation and annihilation 
operators of this free state. So our initial conditions are that h^,^ and its 
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first time derivative coincide with those of at t = 0, and also that 

ipi{x) coincides with ipi{x). The zeroth order solutions to the Heisenberg 
field equations take the form, 

r d^~^k ^ - ^ 

+e;,(r7;^,A)e-*V(fc,A)}, (284) 



C(x) = a-^^^ I j^-yj:-^j:{-^u,ik,s)e^'-bik,s) 



+ -j=Vi{k,X)e~'^-''c^ik,s)} . (285) 

The graviton mode functions are proportional to Hankel functions whose 
precise specification we do not require. The Dirac mode functions Ui{k, s) and 
Vi{k, s) are precisely those of fiat space by virtue of the conformal invariance 
of massless fermions. The canonically normalized creation and annihilation 
operators obey, 

a{k,X),a\k',X')] = 6xy{2nf-^6^-\k-k') , (286) 
{b{k, s), b\k', s')} = 6sA27rf-^6''-\k-k') = {c{k, s), c^k', s')} . (287) 

The zeroth order Fermi field {x) is an anti-commuting operator whereas 
the mode function '^^{x\ k, s) is a (D-number. The latter can be obtained from 
the former by anti-commuting with the fermion creation operator, 

p—ikr) ^ 

v^°(x; fc, s) = a^{^°(x), b\k, s)} = ^u,{k, s)e''-' . (288) 



2k 

The higher order contributions to ipi{x) are no longer linear in the creation 
and annihilation operators, so anti-commuting the full solution ipi{x) with 
b'^lk, s) produces an operator. The quantum-corrected fermion mode function 
we obtain by solving Equation 12561 is the expectation value of this operator 
in the presence of the state which is free vacuum at t = 0, 

^i(x;fc,s) = a^(n\{il}i{x),b^{k,s)}\n) . (289) 

This is what the Schwinger-Keldysh field equations give. The more familiar, 
in-out effective field equations obey a similar relation except that one defines 
the free fields to agree with the full ones in the asymptotic past, and one 
takes the in-out matrix element after anti-commuting. 
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4.3 A Worked-Out Example 

It is perhaps worth seeing a worked-out example, at one loop order, of the 
relation 12891 between the Heisenberg operators and the Schwinger-Keldysh 
field equations. To simplify the analysis we will work with a model of two 
scalars in fiat space, 

L = -d,^*d^v - m W - Ax : V*^ ■■ -^d.xd^x ■ (290) 

In this model (f plays the role of our fermion ipi, and x plays the role of 
the graviton h^^. Note that we have normal-ordered the interaction term 
to avoid the harmless but time-consuming digression that would be required 
to deal with x developing a nonzero expectation value. We shall also omit 
discussion of counterterms. 

The Heisenberg field equations for Equation 12901 are, 

d\~X:^*^: = 0, (291) 
(9^ - m^)<^ - Axv? = . (292) 

As with Dirac + Einstein, we solve these equations perturbatively, 

X{x) = x'{x) + Xx\x)+XWx) + ... , (293) 
ip{x) = ^°(x)+Ay?^(x)+AV(^) + --- • (294) 

The zeroth order solutions are. 

Here k = \\k\\ and u = y/k'^ + m?. The creation and annihilation operators 
are canonically normalized, 

a{k), a^{k') 

We choose to develop perturbation theory so that all the operators and their 
first time derivatives agree with the zeroth order solutions at t = 0. The first 
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b{k),b\k') = c{k),c\k') = (27r)^-^5^-i(^-^') • (297) 



few higher order terms are, 
X\x)= J^dt'J d^-V 

ip\x)=fdti d^'^x'ix 



f}2 / ]-et 
1 



d'^—m? 
1 



d'^ — m? 



(298) 

x') /{x')^\x') , (299) 

/ ret 

x\^J^x\x>\x') + x\x')^\x)]. (300) 



The commutator of ({P{x) with 6"I'(A;) is a C-number, 



V.°(x),6t(^) 



-iujt 



(301) 



However, commuting the full solution with b'^lk) leaves operators. 



^(x), bUk)] = <l>°(x; k) + X f'dt'f d^-^x' (x ^ , x'\ x°(a;')'^°(a;'; fc) 

J /n 7 \ ry^ — 7T),^ /ret 



+>^^ j^dt'l d^'-^x'ix 



d'^ — m? 



X 



ret 



X'(x'),6^(^)]<^°(x') + x'(a;')'^°(a;';^) 
+X°(x')fv^^(x'),6n^)ll+0(A3). (302) 



The commutators in Equation 13021 are easily evaluated. 



x\x'),h\k)y{x') 
t' 



dt" d^'-'x" x' 



1 



x")^^y*{x")ip''ix')^''ix";k) , (303) 



X'{x')[f\x'),b\k) 

rt' 



dt" d^'-'x" x' 



d'^ — m? 



x")^y{x')x\x")^\x"-k).{m 



Hence the expectation value of Equation 13021 gives. 



{p\ [ip{x), h\k)\ \n) = $°(x; k) + \^ j^dt'j d^-^x' (x 

^ ^''"^"{(^1^k">ret(^k°*^^'>°(^') 



d'^ — m? 



X 



ret 



+ (x 



x") (n 

ret 



X"(x')x"(x") n) $"(x"; k) + O(A^). (305) 
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To make contact with the effective field equations we must first recognize 
that the retarded Green's functions can be written in terms of expectation 
values of the free fields, 



X 



Q2 



= -ld{t'-t" 



Q) - (Q 



X 



d'^ — m? 



X 



-ie{t'-t" 



/ / \ 0* / II ^ 



(306) 

X\x")x\x')\n)^ , (307) 
(308) 



-ie{t'-t")\ ln\^\x')^'^*{x")\n) - ln\^*%x")^'^{x')\n) \. (309) 



Substituting these relations into Equation 13051 and canceling some terms gives 
the expression we have been seeking, 



{p\ [p{x), h\k)\ = $°(x; k) - iX^ fdt'l d^-^x' (x 



d'^—m? 



X 



ret 



X'{x')x\x'')\d){^y{x')^'\x'')\d) 



X''{x")x\x') ^'*ix")if'ix') n) $"(x"; k) + 0(A^).(310) 



We turn now to the effective field equations of the Schwinger-Keldysh 
formalism. The C-number field corresponding to ip{x) at linearized order is 
$(2;). If the state is released at t = then the equation $(a;) obeys is, 

(a^ - m2)$(x) - J^dt'J d^-^x'{Ml^{x; x') + Ml_{x; x')}$(x') = . (311) 

The one loop diagram for the self-mass- squared of (p is depicted in Figure |H 
Because the self-mass-squared has two external lines, there are 2^ = 4 
polarities in the Schwinger-Keldysh formalism. The two we require are [151 

El], 



-iMl^{x;x') 



-iXflx 



Q2 



X ) {X 



d'^ — m? 



-iMl_{x-x') = {-i\){+i\)(^x\^\x')^_(^x 



x) + O(A^) , (312) 
+ 0(A^) .(313) 



X 
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\ 



X 



X 



Figure 4: Self-mass-squared for at one loop order. Solid lines stands for 
propagators while dashed lines represent x propagators. 



To recover Equation 13101 we must express the various Schwinger-Keldysh 
propagators in terms of expectation values of the free fields. The ++ polarity 
gives the usual Feynman propagator [81] . 



X 



X 



Q2. 



X = e{t-t'){Q x'{x)x"{x') Q) + e{t'-t){Q x'{x')x'{x) ,(314) 



X 



= e{t-t')(Q\<f\x)if'^*{x')\Q)+e{t'-t)(Q\if^*{x')^\x)\Q). (315) 
The +- polarity propagators are [81] . 



X 



Q2 



X 



X 



x\x)x\x) n) , 



^'*{x')v\x) 



(316) 
(317) 



Substituting these relations into Equation 13121 and Equation 13131 and making 
use of the identity 1 = 6{t — t')+6{t'—t) gives, 



Ml^{x;x') + Ml_{x;x') = -tX^e{t-t')\ (r]|x°(a;)x°(a;') 



x(n 



y^°(x)(/^°*(x')|^^)-(^^|x°(x')x°(a;)|^^)(^^|</'°*(a;')¥^°(x)|fi)l+0(A^).(31 
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We now solve Equation 13111 perturbatively. The free plane wave mode 
function |3U II is of course a solution at order A°. With Equation 13181 we easily 
recognize its perturbative development as, 



l^dt'J d^-^x' (x 


1 







X 



ret 



-in 



X\x")x\x') 



n) |>$"(x"; k) + 0(A^).(319) 



That agrees with Equation I310[ so we have established the desired connec- 
tion, 



A;) = (n ^{x),b\k) 



(320) 



at one loop order. 



4.4 Gauge Issues 

The preceding discussion has made clear that we are working in a particular 
local Lorentz and general coordinate gauge. We are also doing perturbation 
theory. The function k, s) describes how a free fermion of wave number 
k and helicity s propagates through classical de Sitter background in our 
gauge. What s) gives is the first quantum correction to this mode 

function. It is natural to wonder how the effective field \l/j(a;; k, s) changes if 
a different gauge is used. 

The operators of the original, invariant Lagrangian transform as follows 
under diffeomorphisms {x^ — > x'^) and local Lorentz rotations (Ajj)|f| 

^:(a;) = A,,{x'-\x))^P,{x'-\x)) , (321) 

eU(x) = ^A,^{x'-\x))e^^{x'-\x)) . (322) 

^Of course the spinor and vector representations of the local Lorentz transformation 
are related as usual, with same parameters LUcd{x) contracted into the appropriate repre- 
sentation matrices, 

Ay EE 5^j - ^ujcdJ% + ■■■ and A^^^ ^ S^"" - uj^" + . . . . 
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The invariance of the theory guarantees that the transformation of any solu- 
tion is also a solution. Hence the possibility of performing local transforma- 
tions precludes the existence of a unique initial value solution. This is why no 
Hamiltonian formalism is possible until the gauge has been fixed sufficiently 
to eliminate transformations which leave the initial value surface unaffected. 

Different gauges can be reached using field-dependent gauge transforma- 
tions [82j. This has a relatively simple effect upon the Heisenberg operator 
ipi{x), but a complicated one on the linearized effective field \l/i(a;; k, s). Be- 
cause local Lorentz and diffeomorphism gauge conditions are typically spec- 
ified in terms of the gravitational fields, we assume x''^ and Aj^ depend upon 
the graviton field h^^. Hence so too does the transformed field, 

ij'Mix) = KAh]{x'~'[h]{x))4j,{x'~'[h]{x)) . (323) 

In the general case that the gauge changes even on the initial value surface, 
the creation and annihilation operators also transform, 

b'[h]{k,s) = -^^:ik,s) J d^'-'xe-^'-^Mr^^x) , (324) 



where 7]i = —1/H is the initial conformal time. Hence the linearized effective 
field transforms to, 

'^■{x;k,s) =a^(^Q\[^[[h]{x),b'^[h]{k,s)}\Q) . (325) 

This is quite a complicated relation. Note in particular that the /i^j^ de- 
pendence of x'^[h] and Ajjj/i] means that \l/'j(x; k, s) is not simply a Lorentz 
transformation of the original function \l/j(x;A;,s) evaluated at some trans- 
formed point. 



5 ENHANCED FERMION MODE FUNC- 
TION 

We first modify our regularized result for the fermion self energy by the 
employing Schwinger-Keldysh formalism to make it causal and real. We 
then solve the quantum corrected Dirac equation and find the fermion mode 
function at late times. Our result is that it grows without bound as if there 
were a time-dependent field strength renormalization of the free field mode 
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function. If inflation lasts long enough, perturbation theory must break 
down. The same result occurs in the Hartree approximation although the 
numerical coefficients differ. 



5.1 Some Key Reductions 

The purpose of this section is to derive three results that are used repeatedly 
in reducing the nonlocal contributions to the effective field equations. We 
observe that the nonlocal terms of Equation 12511 contain l/Ax^. We can 
avoid denominators by extracting another derivative, 

— — = — In Ax^ and , , ' = — In^ Ax^ - 2 In Ax' 

Ax2 4 ^ ^ Ax2 8 L ^ ^ ^ 

(326) 

The Schwinger-Keldysh field equations involve the difference of ++ and +- 
terms, for example, 

Ax2^ Ax2_ 
= y|ln2(/i2Ax2^) -21n(/i2Ax^^) -ln2(/i2Ax2j + 21n(/i2Ax2j|.(327) 

We now define the coordinate intervals Arj = rj — r]' and Ax = ||x — x'|| in 
terms of which the ++ and +- intervals are, 

Ax^^ = Ax^ - {\Ar]\-i6Y and Ax^_ = Ax^ - {Ar]+i6y . (328) 

When 1]' > rj we have Ax^_|_ = Ax^_ , so the ++ and +- terms in Equation 13271 
cancel. This means there is no contribution from the future. When rj' < r) 
and Ax > At] (past spacelike separation) we can take 6 = 0, 

ln{fj.^Axl^) = ln[/i2(Ax2-Ar/2)] = ln{fi^Axl_) {Ax > At] > 0) . (329) 

So the ++ and +- terms again cancel. Only for rj' < rj and Ax < Arj (past 
timelike separation) are the two logarithms different, 

\n{i?Axl^) = ln[fi^{Arj^-Ax^)] ± in (Arj > Ax > 0) . (330) 

Hence Equation 13271 can be written as, 

ln(^2^x^^^) _ Hf^^Axl_)_ ^ ^JLd^UArj-Ax)\HA^v'-^x-)-l]]-im 



Ax^ Ax? 
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This step shows how the Schwinger-Kledysh formahsm achieves causahty. 

To integrate expression 133 1 1 up against the plane wave mode function 12811 
we first pull the x'^ derivatives outside the integration, then make the change 
of variables x' = x+r and perform the angular integrals, 



A ln(/i2Ax^+) \xi{^'^Ax\ 



vl/°(V,x',A;,s) 



Axi 



•J Tli \/ An, ^ 



2k Jo 

V 



kV2k 



Vi 



X 



j\izz sin(a2){ln(l -^2) + 2 In(^) - 1}. (332) 



Here a = kAr] and r]i = — is the initial conformal time, corresponding to 
physical time t = 0. The integral over z is facilitated by the special function, 



1 2 1 
dzz sm{az) ln(l — 2;^) = — sin(a) - 



cosfaj + asinfaj 



si(2a) + | 



1 r 
+ — 



sm a — acos a 



ci(2a)-7-ln(-)J. (333) 



Here 7 is the Euler-Mascheroni constant and the sine and cosine integrals 
are, 



SI X 



ci X = 



dt 



sin(t) 



TT 



t 



-+ dt 



sint 







t 



cost , , , , , 

dt — [ — = 7 + ln(x) + dt 



t 







cos(t) — 1 



t 



(334) 
(335) 



After substituting the ^ function and performing the elementary integrals. 
Equation 13321 becomes. 



d X 



Xj 



ln(/i^Ax^ 



Axi 



mv',x,k,s) 



z27r2 



Ak-x 



kV2k 



Ui{k, s) 



xidi^ + l)ldv'e-^'^{a'ac^)+ 21n( 



-1 



sm a —a cos a 



-. (336) 



One can see that the integrand is of order a^ln(a) for small a, which 
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means we can pass the derivatives through the integraL After some rear- 
rangements, the first key identity emerges, 



i47r^/c"^^'"(?7, x; k, s) / drj'e 



J ik Arj 



cos{kAri) / dt — — 



+ sin(A;Ar7) 



2fcA^ COs(t)-l 
V t 



-21n(2/iA?7) 



(337) 



Note that we have written e = e x e^^^^"^ and extracted the first phase 
to reconstruct the full tree order solution '^^{rj, x; k, s) = ^^-^Ui{k, s)e**^'^. 

The second identity derives from acting a d'Alembertian on Equation [3371 
The d'Alembertian passes through the tree order solution to give. 



d'^^\r],x; k,s) = ^\'r],x;k, s)d^{dr,-2ik) . 



(338) 



Because the integrand goes like aln(a) for small a, we can pass the first 
derivative through the integral to give. 



/ d^x 



,j\n{ij^Axl^) ln(/i2Ax^ 



^°(r/',f' ; A;,s) 



Axl 



= «47r'*°(r7,f;^,s)9^j^V|^'"t^^(^^) +21n(^)|. (339) 

We can pass the final derivative through the first integral but, for the second, 
we must carry out the integration. The result is our second key identity, 

J Ax^, , Axi 



i47r^^"(r/,f;fc,s)<^21n 



(340) 



The final key identity is derived through the same procedures. Because 
they should be familiar by now we simply give the result. 



1 



1 



^°(V,f ;fc,s) 



Axl^ 



iAn'k~^'^%ri,x;k,s) I dr]' e'"'^'^ sm{kAri) . (341) 



/ ikArj , 



Vi 
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Table 19: Derivative operators f//-: Their common prefactor is ^ 



/ 




I 




1 




4 


-8 


2 




5 


4 ^V^ 


3 




6 


7 ^V' 



5.2 Solving the Effective Dirac Equation 



In this section we first evaluate the various nonlocal contributions using the 
three identities of the previous section. Then we evaluate the vastly simpler 
and, as it turns out, more important, local contributions. Finally, we solve 
for '^^{t], X] k, s) at late times. 

The various nonlocal contributions to Equation 12561 take the form, 



1=1 



J j ln(af A 



\n{ajAxl ) 1 T-o/ / -/ r n 
A^2 ^-^ j^',iv',x'-k,s) 



+ id^x'Uf, 



1 



Ax2 



.vl/0(r^',f';fc,s) . (342) 



The spinor differential operators Ul^ are listed in Table [191 The constants 
are /i for J = 1, 2, 3, and \H for / = 4, 5. 



As an example, consider the contribution from Uf,: 



15 K^H^ 
T 2^4 



,|ln(/i2Aa;^^) ^(/i^Aa;^ 



Ax2^ 



15 K^H^ 



2 2%4 



^xz47r2^°(77,f;fc, s)|21n 
15 1 



Ax2 
-2^ 



l + /7r/) +/ dr]' 



iHY'^ iv,^', k, s) X 



2 1+ifr/ 



,2i-|(l+//r,) 



>,(343) 
(344) 



In these reductions we have used i ^\E'°(r7, x; k, s) = i7°\E'°(r7, x; k, s) and 
the second key identity 13401 Recall from the Introduction that reliable pre- 
dictions are only possible for late times, which corresponds to r/ ^ 0^. We 
therefore take this limit, 



15 K^H^ 
2 2^4 



d^x- 



,|ln(/i2Ax^+) ln(/i2Ax^_ 



Axl^ 



Axl_ 
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.^°(V,f ;fc,s) 



Table 20: Nonlocal contributions to J d^x'[T]{x;x')'^^{ri',x';k,s) at late 
times. Multiply each term by x iH'j'^'^^^r], x; k,s). 



T 


r^opfln PlP'^"^ of "f"Vip Ipfp "f"iTnp PO'^"^T'll'~nl"^^o'n fyrtTn p^ipTi TI^ 

'J 


1 





2 


f {exp(2z|) + l} 


3 


-z|{21n(t)-/°^r^-(-P(^y-'')-^)} 


4 




5 


4f /° rfVe-^^'='''{/o-''''' + 2 ln(ifr^')} 


6 


-|z|{exp(2z|) - 1} 



— > ^^H^^^>\r^,x■k,s) x ^{exp(2z-|) + l}.(345) 

The other five nonlocal terms have very similar reductions. Each of them 
also goes to ^i^^ x zif7°\l/°(?7, x; k, s) times a finite constant at late times. We 
summarize the results in Table [20] and relegate the details to an appendix. 

The next step is to evaluate the local contributions. This is a straightfor- 
ward exercise in calculus, using only the properties of the tree order solution 
12811 and the fact that (9^a = Ha?5^^. The result is, 

,2 tt2 



J y H^aa' 2 J 



+ y (ln(a) ^ ^ln(a))^°(r/, x; fc, s) - 14 ln(a) ^^^(r], f; k, (346) 
t^^H"^ ... n.n. ^ . fir . A; , _ „. " 



2g^2^ff7"^"(^,^;/^,s) X |_a-14z-ln(a)-2^-|. (347) 

The local quantum corrections 13471 are evidently much stronger than their 
nonlocal counterparts in Table [201 Whereas the nonlocal terms approach 
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a constant, the leading local contribution grows like the inflationary scale 
factor, a = e^*. Even factors of ln(a) are negligible by comparison. We can 
therefore write the late time limit of the one loop fleld equation as, 

i^K''-^\r],x;k,s) ^iJL}l,Ha-f''^I/\r],x;k,s) . (348) 

The only way for the left hand side to reproduce such rapid growth is for the 
time derivative in to act on a factor of ln(a), 

i^d^. Ha) = iY^^l = iHa-f'^ . (349) 

We can therefore write the late time limit of the tree plus one loop mode 
functions as, 

^ f K^H^ 17 1 

'$^{r],x;k,s) + K'^'^\r],x;k,s) — > \ 1 + — — — ln(a) \^°{r],x;k, s) . (350) 

[ /"vr^ / J 

All other corrections actually fall off at late times. For example, those from 
the ln(a) terms in Equation 13471 go like ln(a)/a. 

There is a clear physical interpretation for the sort of solution we see 
in Equation 13501 When the corrected fleld goes to the free fleld times a 
constant, that constant represents a fleld strength renormalization. When 
the quantum corrected fleld goes to the free fleld times a function of time 
that is independent of the form of the free fleld solution, it is natural to think 
in terms of a time dependent field strength renormalization, 

vl>(r/,f;^,s)-^^^^^^^lflM where Z^it) = l-^^Ha) + 0{K') . 

(351) 

Of course we only have the order correction, so one does not know if 
this behavior persists at higher orders. If no higher loop correction su- 
pervenes, the fleld would switch from positive norm to negative norm at 
ln(a) = 2^7r^/17K^i/^. In any case, it is safe to conclude that perturbation 
theory must break down near this time. 

5.3 Hartree Approximation 

The appearance of a time- dependent fleld strength renormalization is such 
a surprising result that it is worth noting we can understand it on a simple. 
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qualitative level using the Hartree, or mean-field, approximation. This tech- 
nique has proved useful in a wide variety of problems from atomic physics 
[83] and statistical mechanics [84], to nuclear physics [85] and quantum field 
theory [86j. Of particular relevance to our work is the insight the Hartree 
approximation provides into the generation of photon mass by inflationary 
particle production in SQED [STllHSlEn]. 

The idea is that we can approximate the dynamics of Fermi fields inter- 
acting with the graviton field operator, h^u, by taking the expectation value 
of the Dirac Lagrangian in the graviton vacuum. To the order we shall need 
it, the Dirac Lagrangian is Equation [50l 



+K^\-h'^--hP''hp, 



Of course the expectation value of a single graviton field is zero, but the 
expectation value of the product of two fields is the graviton propagator in 
Equation [711 



{Vt I T h^y{x)hpa{x) I Vt) 



i/S.A(x\ x') 



rpA 

t^l^''' per 



+ iAsix; x') 



+ iAc{x; x') 



rnC 
fJ.'^-'' pa 



(353) 



Recall the index factors from Equations [76]|78[ 

2 



lif-'- pa 



'^■n^.i^pTl^y 



D-3 



VfJ.i'Vpa 



^f^-^ pa 



) ' 



rpC 

^^-^ pa 



(D-2)(D-3) 



(D-3)« + r/^, {D-?,)6X + npa 



(354) 
(355) 



Recall also that parenthesized indices are symmetrized and that a bar over 
a common tensor such as the Kronecker delta function denotes that its tem- 
poral components have been nulled. 



5^. = 5^. - (5,^5.^ 



Vpu = Vp. + ■ (356) 



The three scalar propagators that appear in Equation 13531 have compli- 
cated expressions 18311851 which imply the following results for their coincidence 
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limits and for the coincidence limits of their first derivatives, 



lim iAA^x; x') 
lim d^iAA^x; x) 

a;'— >x 

lim iAb^x; x) 
lim d^iAsix; x') 

a;'— >x 

lim iAc{x; x) 
lim d^j_iAc{x]x') 



H 



D-2 



r(D-i) 



2 



r(f; 



HD-2 y(D - 1] 

(47r)f r(f) ' ^ i'^ 

HD-2 y(D - 1) 1 



-7rcot(|D) +21n(a) !>, (357) 



X HaS^ = lim d'iAAix; x') , (358) 



(47r)f r(f; 
T{D-r 



X 



D-2 



0= lim d'^iAsixix') 



H 



x'-^x 
D-2 



X 



(47r)f r(f) - {D-2){D-3) 
= lim d'^iAc{x; x) . 



(359) 
(360) 

(361) 
(362) 



We are interested in terms which grow at late times. Because the 5-type and 
C-type propagators go to constants, and their derivatives vanish, they can 
be neglected. The same is true for the divergent constant in the coincidence 
limit of the A-type propagator. In the full theory it would be absorbed into 
a constant counterterm. Because the remaining, time dependent terms are 
finite, we may as well take Z) = 4. Our Hartree approximation therefore 
amounts to making the following replacements in Equation 1352^ 



IP 

8^^ 



/la li/p 



VppVua + V^aVup-^Vf^uV, 



pa I up 



Ipu I pcj 



(363) 
(364) 



It is now just a matter of contracting Equations 1363113641 appropriately to 
produce each of the quadratic terms in Equation 13521 For example, the first 
term gives. 



2^2 
2^2 



ln(a) 
ln(a) 



3 + 3-1 



VppVua + ripaVyp - 2rip^ripJ\^-<^i^'<^, (365) 
Wi^^ . (366) 



The second quadratic term gives a proportional result, 

-2 -k'H^ 



— K 



2%2 



ln(a)[9 + 3-6 ^i^^ . 



(367) 
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The total for these first two terms is ~^^2^^ ln(a)\l/i^\l/. 

The third and fourth of the quadratic terms in Equation 13521 result in 
only spatial derivatives, 



— K 



- ln(a' 
ln(a) 



1 + 1-6 
3 + 1-2 



(368) 
(369) 



The total for this type of contribution is ln(a)\I'i(p\E'. 

The final four quadratic terms in Equation 13521 involve derivatives acting 
on at least one of the two graviton fields, 



4 



K 



-Ha 



K 



2%2 
. 



Ha 
Ha 



1 + 1-6 
3 + 1-2' 
3 + 1-2' 



r7^p^7'V^°^ , (370) 

77^,^7°^''"^, (371) 
r/^,^7^J^°^, (372) 
(373) 



2 ' 



The second of these contributions vanishes owing to the antisymmetry of the 
Lorentz representation matrices, J^" = |[7'^,7'^], whereas ^/^p7^J''° ■ 
Hence the sum of all four terms is Ha^i'y'^^. 
Combining these results gives. 



Dirac 



— 3^2^2 _ 

^^^^ ln(a)\l/z^\Ef 



4vr2 



-Ha'^iY'^ + 



3^2 
87r2 



■ ln(a) 



3k^H^ 



ln(a) 



^4 



167r2 
7k^H^ 



ln(a)^i^^ + 0(K*), (374) 



167r2 



ln(a)*2^^+0(K^).(375) 



If we express the equations associated with Equation 13751 according to the 
perturbative scheme of Section 2, the first order equation is, 

i^K^^\T],x]k,s) = —-^iH-f°'^\r],x]k,s){2Aa-28i—ln{a)} . (376) 
2 vr ^ rf J 



This is similar, but not identical to, what we got in expression 13471 from 
the delta function terms of the actual one loop self-energy in Equation 12511 
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In particular, the exact calculation gives y'^" 14?-|: ln(a), rather than the 
Hartree approximation of 24a — 28i-|: ln(a). Of course the ln(a) terms make 
corrections to \E'^ which fall like ln(a)/a, so the real disagreement between 
the two methods is limited to the differing factors of ^ versus 24. 

We are pleased that such a simple technique comes so close to recovering 
the result of a long and tedious calculation. The slight discrepancy is no 
doubt due to terms in the Dirac Lagrangian by Equation 13521 which are 
linear in the graviton field operator. As described in relation [2891 of section 2, 
the linearized effective field \E'i(x; k, s) represents times the expectation 
value of the anti-commutator of the Heisenberg field operator ipi{x) with the 
free fermion creation operator b{k, s). At the order we are working, quantum 
corrections to \E'j(x; k, s) derive from perturbative corrections to ipii^) which 
are quadratic in the free graviton creation and annihilation operators. Some 
of these corrections come from a single hhipip vertex, while others derive 
from two hipip vertices. The Hartree approximation recovers corrections of 
the first kind, but not the second, which is why we believe it fails to agree 
with the exact result. Yukawa theory presents a fully worked-out example 
pTl [T2l [90] in which the entire lowest-order correction to the fermion mode 
functions derives from the product of two such linear terms, so the Hartree 
approximation fails completely in that case. 

6 CONCLUSIONS 

We have used dimensional regularization to compute quantum gravitational 
corrections to the fermion self-energy at one loop order in a locally de Sit- 
ter background. Our regulated result is Equation 12381 Although Dirac + 
Einstein is not perturbatively renormalizable [18j we obtained a finite result 
shown by Equation 1251! by absorbing the divergences with BPHZ countert- 
erms. 

For this IPI function, and at one loop order, only three counterterms are 
necessary. None of them represents redefinitions of terms in the Lagrangian 
of Dirac + Einstein. Two of the required counterterms of Equation [^2] are 
generally coordinate invariant fermion bilinears of dimension six. The third 
counterterm of Equation 1104! is the only other fermion bilinear of dimension 
six which respects the symmetries shown by Equations [6T][66[ of our de Sit- 
ter noninvariant gauge shown in Equation [SD] and also obeys the reflection 
property shown in Equation 11031 of the self-energy for massless fermions. 
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Although parts of this computation are quite intricate we have good con- 
fidence that Equation 12511 is correct for three reasons. First, there is the 
fiat space hmit of taking H to zero while taking the conformal time to be 
7] = —e~^^/H with t held fixed. This checks the leading conformal contribu- 
tions. Our second reason for confidence is the fact that all divergences can 
be absorbed using just the three counterterms we have inferred in chapter 
2 on the basis of symmetry. This was by no means the case for individual 
terms; many separate pieces must be added to ehminate other divergences. 
The final check comes from the fact that the self-energy of a massless fermion 
must be odd under interchange of its two coordinates. This was again not 
true for separate contributions, yet it emerged when terms were summed. 

Although our renormalized result could be changed by altering the finite 
parts of the three BPHZ counterterms, this does not affect its leading behav- 
ior in the far infrared. It is simple to be quantitative about this. Were we to 
make finite shifts Aa^ in our counterterms Equation 12471 the induced change 
in the renormalized self-energy would be. 



—I 



[x: x') = —K^< 



^'^-^^d^ + l2Aa2H^^+Aa:iH^^\5\x--x') . (377) 



aa 



No physical principle seems to fix the Aoj so any result that derives from 
their values is arbitrary. This is why BPHZ renormalization does not yield 
a complete theory. However, at late times (which accesses the far infrared 
because all momenta are redshifted by a{t) = e^*) the local part of the renor- 
malized self-energy of Equation 12511 is dominated by the large logarithms, 

[In(aa') + l^i^(^aa')H^^ - 7\n{aa')H'A6\x-x') . (378) 



267r2 [ aa' 

The coefficients of these logarithms are finite and completely fixed by our 
calculation. As long as the shifts Aa, are finite, their impact Equation 13771 
must eventually be dwarfed by the large logarithms in Equation 13781 

None of this should seem surprising, although it does with disturbing 
regularity. The comparison we have just made is a standard feature of low 
energy effective field theory and has a very old and distinguished pedigree 
[2aiM[2a[26l[27l[28l[29llMll3ll|33[3^ Loops of massless 

particles make finite, nonanalytic contributions which cannot be changed by 
local counterterms and which dominate the far infrared. Further, these effects 
must occur as well, with precisely the same numerical values, in whatever 
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fundamental theory ultimately resolves the ultraviolet problem of quantum 
gravity. That is why Feinberg and Sucher got exactly the same long range 
force from the exchange of massless neutrinos using Fermi theory [251 126] as 
one would get from the Standard Model [26] . 

So we can use Equation 12511 reliably in the far infrared. Our motivation 
for undertaking this exercise was to search for a gravitational analogue of 
what Yukawa-coupling a massless, minimally coupled scalar does to massless 
fermions during inflation [11]. Obtaining Equation 12511 completes the first 
part in that program. In the second stage we used the Schwinger-Keldysh 
formalism to include one loop, quantum gravitational corrections to the Dirac 
equation. Because Dirac + Einstein is not perturbatively renormalizable, 
it makes no sense to solve this equation generally. However, the equation 
should give reliable predictions at late times when the arbitrary finite parts 
of the BPHZ counterterms Equation 12461 are insignificant compared to the 
completely determined factors of In(aa') on terms of Equations I258II260I which 
otherwise have the same structure. In this late time limit we find that the 
one loop corrected, spatial plane wave mode functions behave as if the tree 
order mode functions were simply subject to a time-dependent field strength 
renormalization, 

Z^it) = 1 - —GHHn(a) + 0(G^) where G = Wttk^ . (379) 
An 

If unchecked by higher loop effects, this would vanish at ln(a) ~ l/GH"^. 
What actually happens depends upon higher order corrections, but there is 
no way to avoid perturbation theory breaking down at this time, at least in 
this gauge. 

Might this result be a gauge artifact? One reaches different gauges by 
making field dependent transformations of the Heisenberg operators. We 
have worked out the change in Equation 13251 this induces in the linearized 
effective field, but the result is not simple. Although the hnearized effective 
field obviously changes when different gauge conditions are employed to com- 
pute it, we believe (but have not proven) that the late time factors of ln(a) 
do not change. 

It is important to realize that the IPI functions of a gauge theory in a 
fixed gauge are not devoid of physical content by virtue of depending upon 
the gauge. In fact, they encapsulate the physics of a quantum gauge field 
every bit as completely as they do when no gauge symmetry is present. One 
extracts this physics by forming the IPI functions into gauge independent 
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and physically meaningful combinations. The S-matrix accomplishes this in 
flat space quantum fleld theory. Unfortunately, the S-matrix fails to exist 
for Dirac + Einstein in de Sitter background, nor would it correspond to an 
experiment that could be performed if it did exist [HI |92l 193] . 

If it is conceded that we know what it means to release the universe in a 
free state then it would be simple enough — albeit tedious — to construct an 
analogue of ipilx) which is invariant under gauge transformations that do not 
affect the initial value surface. For example, one might extend to fermions 
the treatment given for pure gravity by [9^ : 

• Propagate an operator-valued geodesic a fixed invariant time from the 
initial value surface; 

• Use the spin connection A^cdJ'^'^ to parallel transport along the geodes- 
ic; and 

• Evaluate if) at the operator-valued geodesic, in the Lorentz frame which 
is transported from the initial value surface. 

This would make an invariant, as would any number of other constructions 
[95j . For that matter, the gauge-fixed IPI functions also correspond to the 
expectation values of invariant operators [82]. Mere invariance does not guar- 
antee physical significance, nor does gauge dependence preclude it. 

What is needed is for the community to agree upon a relatively simple 
set of operators which stand for experiments that could be performed in de 
Sitter space. There is every reason to expect a successful outcome because 
the last few years have witnessed a resolution of the similar issue of how to 
measure quantum gravitational back-reaction during infiation, driven either 
by a scalar infiaton [HSl EZl ESI ES] or by a bare cosmological constant [100] . 
That process has begun for quantum field theory in de Sitter space [HIl [121 
|95| 1100] and one must wait for it to run its course. In the meantime, it is 
safest to stick with what we have actually shown: perturbation theory must 
break down for Dirac + Einstein in the simplest gauge. 

This is a surprising result but we were able to understand it qualitatively 
using the Hartree approximation in which one takes the expectation value of 
the Dirac Lagrangian in the graviton vacuum. The physical interpretation 
seems to be that fermions propagate through an effective geometry whose 
ever-increasing deviation from de Sitter is controlled by infiationary graviton 
production. At one loop order the fermions are passive spectators to this 
effective geometry. 
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It is significant that inflationary graviton production enhances fermion 
mode functions by a factor of ln(a) at one loop. Similar factors of ln(a) have 
been found in the graviton vacuum energy [651 166]. These infrared logarithms 
also occur in the vacuum energy and mode functions of a massless, minimally 
coupled scalar with a quartic self- interaction |56l |571 1101] , and in the VEV's 
of almost all operators in Yukawa theory [20] and SQED [W2\ 1103] . A recent 
all orders analysis was not even able to exclude the possibility that they 
might contaminate the power spectrum of primordial density fluctuations 
[TM[T05lfT06] ! 

The fact that infrared logarithms grow without bound raises the exciting 
possibility that quantum gravitational corrections may be significant dur- 
ing inflation, in spite of the minuscule coupling constant of GH"^ < 10~^^. 
However, the only thing one can legitimately conclude from the perturbative 
analysis is that infrared logarithms cause perturbation theory to break down, 
in our gauge, if inflation lasts long enough. Inferring what happens after this 
breakdown requires a nonperturbative technique. 

Starobinskii has long advocated that a simple stochastic formulation of 
scalar potential models serves to reproduce the leading infrared logarithms 
of these models at each order in perturbation theory [107] . This fact has 
recently been proved to all orders [1081 1109j . When the scalar potential is 
bounded below it is even possible to sum the series of leading infrared loga- 
rithms and infer their net effect at asymptotically late times [110 ] ! Applying 
Starobinskii's technique to more complicated theories which also show in- 
frared logarithms is a formidable problem, but solutions have recently been 
obtained for Yukawa theory [90] and for SQED [103] . It would be very inter- 
esting to see what this technique gives for the infrared logarithms we have 
exhibited, to lowest order, in Dirac + Einstein. And it should be noted 
that even the potentially complicated, invariant operators which might be 
required to settle the gauge issue would be straightforward to compute in 
such a stochastic formulation. 



7 NONLOCAL TERMS FROM TABLE 5.2 



It is important to establish that the nonlocal terms make no significant con- 
tribution at late times, so we will derive the results summarized in Table [201 
For simplicity we denote as [U^] the contribution from each operator U^j in 

Table [T9l We also abbreviate \l/'^(r7,a;; k,s) as \E''^(x). 
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Owing to the factor oil /a' in f/?^-, and to the larger number of derivatives, 
the reduction of \U^] is atypical, 



1 



—IK 



/ (fx 



, 1 [ ln(/i2Ax^+) ln(/i2Ax^ 



1 Axt 



Ax2 



7°*°(x) [-2zfc9, + d^^ I 9, rd7]'{-H7]') ( 



(380) 



—IK 



+dl / rfV(-2/7r/')ln(2/iAr7) |>, (381) 



K^ 



2 



+ 



(1 + /7r/) (1 + Hr]y 



(382) 
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Hrj 



(383) 



This expression actually vanishes in the late time limit of 77 . 
[[/2] was reduced in Section 4 so we continue with [[/^], 
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,J ln(yu2Ax 
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r2/i 
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(384) 
,(385) 
(386) 
(387) 



Ufj has the same derivative structure as Uf^, so [f/^] follows from Equation 
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K^ 

' 2H^ 
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1 Tj2 

Aa;2^ 



Ax2 



^'\x') , (388) 
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267r2 



8A;7°^°(x)<j21n[(l + ifr7) 



d^'{ )^(389) 



(390) 



C/f- has a Laplacian rather than a d'Alembertian so we use identity 13371 



for \U^]. We also employ the abbreviation kAri = a, 



(391) 



++ 



f , , /-^a, sin()f:) . , ,r/"^", /cos(t)-l 
x< — cos(a;)y at — ^ + sin(a; 



t 



21n( 



t72 iL2 „2a 



-it 



1^2 0-2 i,2 ^0 r /•2q 



-it 



t\2 



+ ln(i7r/') 



(392) 
(393) 
(394) 



Uf, has the same derivative structure as Uf, but it acts on a different inte- 
grand. We therefore apply identity 13411 for [f/ ], 



(395) 



X (-i47r2)A;-i^°(x) / ' rfr^V'^^" sin(fcAr7) , (396) 
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